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■ We describe the effective lagrangian approach to the color superconductivity. 

Superconductivity in color interactions is due to vacuum instability at T = 
Q_i| and high densities, which arises because the formation of quark-quark pairs is 

(~) • energetically favored. The true vacuum is therefore populated by a condensate 

of Cooper pairs breaking color and baryonic number. The effective langrangian 
follows from the Wilson's renormalization group approach and is based on the 
idea that at T = and \x — ► oo the only relevant degrees of freedom are those 
near the Fermi surface. The effective description that arises if one considers 
only the leading terms in the 1/ ' \i expansion is particularly simple. It is based 
on a lagrangian whose effective fermion fields are velocity-dependent; moreover 
strong interactions do not change quark velocity (Fermi velocity supers election 
rule) and the effective lagrangian does not contain spin matrices. All these 
features render the effective theory similar to the Heavy Quark Effective Theory, 
which is the limit of Quantum ChromoDynamics for mg — ► oo. For this reason 
one can refer to the effective lagrangian at high density as the High Density 
Effective Theory (HDET). In some cases HDET results in analytical, though 
approximate, relations that are particularly simple to handle. 

After a pedagogical introduction, several topics are considered. They in- 
clude the treatment of the Color-Flavor-Locking and the 2SC model, with eval- 
uation of the gap parameters by the Nambu-Gorkov equations, approximate 
dispersion laws for the gluons and calculations of the Nambu-Goldstone Bosons 
properties. We also discuss the effective lagrangian for the crystalline color 
superconductive (LOFF) phase and we give a description of the phonon field 
related to the breaking of the rotational and translational invariance. Finally a 
few astrophysical applications of color superconductivity are discussed. 
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Chapter 1 

A pedagogical introduction 



1.1 Overview 

The aim of the present review is to describe various approaches to the color 
superconductivity (CSC) that are based on the method of the effective la- 
grangians. CSC is one of the most fascinating advances in Quantum Chromo 
Dynamics (QCD) in recent years |, |, |, §, g |, |, |, |0|, [jl], 0, [l3j (for 



re- 



views see [|14|, [l5|, |l6|); it offers a clue to the behavior of strong interactions 
at very high baryonic densities, an issue of paramount relevance both for the 
understanding of heavy ion collisions and the physics of compact stars. 

The origin of the superconductivity in color interactions is related to the 
instability of the naive vacuum at T = 0; at low densities this instability is 
supposed to arise because the formation of a quark-antiquark pair is energeti- 
cally favored; therefore Cooper pairs formed by a quark and an antiquark (in 
a color singlet representation) populate the true vacuum. This phase is indeed 
characterized by the order parameter 

(^)^O, (1.1) 

i.e. by a condensate of qq pairs, similarly to the phenomenon of the Cooper 
pairs in ordinary superconducting materials at low temperatures, described by 
the theory of Bardeen, Cooper and Schrieffer (BCS)[17]. The condensate ([O]) 



violates chiral invariance, which is a symmetry of the QCD Lagrangian if one 
neglects the current quark masses. Therefore a non vanishing condensate ( |l.l| ) 
is tantamount to the spontaneous breaking of the chiral symmetry and the 
existence, because of the Goldstone's theorem, of a certain number of massless 
particles, the Nambu-Goldstone-Bosons (NGBs). 

At high baryonic densities a different phenomenon sets in. As a matter 
of fact, for sufficiently high baryon chemical potential fi, the color interaction 
favors the formation of a quark-quark condensate in the color antisymmetric 
channel. This condensate 

(C^/j), (1.2) 

(a, f3 = 1,2,3 color indices; i, j = 1,2,3 flavor indices) acts as the order pa- 
rameter of a new phase where both the SU(3) C color symmetry and the U(1)b 
baryonic number are spontaneously broken. As to chiral symmetry, the color 
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interaction favors in general pairing of quarks of total spin and opposite mo- 
menta; therefore the quarks in (|1.2| ) have the same helicities and the condensate 
(1.2) does not break chiral invariance explicitly^. 

Different phenomena take place depending on the value of the order param- 



eter (1.2). One could have indeed: 

e^-'dj < ifaCfyf) >= A^ 3 (1.3) 

where the sum over the flavor indices run from 1 to 2 and tp represents a left 
handed 2-component Weyl spinor (the right handed field satisfies the same 
relation with A — ► —A); moreover a sum over spinor indices is understood and 
C = ia 2 - This case correspond to the decoupling of the strange quark (m s — > 
oo; m u = rrid = 0) and is called in the literature 2SC model. From dynamical 
analyses [0], |[4] one knows that, for fj, sufficiently large, the condensate ( |1.3| ) 
is non vanishing. Therefore it breaks the original symmetry group SU(3) C <8> 
SU{2) L <g> SU{2) R ® U(1) B down to 

SU(2) C ® SU(2) L <g> SU(2) R g> Z 2 . (1.4) 

The chiral group remains unbroken, while the original color symmetry group is 
broken to SU(2) C , with generators T A corresponding to the generators T 1 , T 2 , T 3 
of SU (3) c . As a consequence, three gluons remain massless whereas the remain- 
ing five acquire a mass. The Z 2 group means that the quark fields can still be 
multiplied by -1. Even though the original U(1)b is broken there is an unbro- 
ken global symmetry that plays the role of U(1)b- As for U(1)a, this axial 
symmetry is broken by anomalies, so that in principle there is no Goldstone 
boson associated to its breaking by the condensate; however at high densities 
explicit axial symmetry breaking is weak and therefore there is a light would 
be Goldstone boson associated to the breaking of the axial color. One can con- 
struct an effective theory to describe the emergence of the unbroken subgroup 
SU(2) C and the low energy excitations, much in the same way as one builds 
up chiral effective lagrangian with effective fields at zero density. For the two 
flavor case this development can be found in Jl8|]. 

For the three flavor case (m u = md = m s = 0) the following interesting 
caseQ has been widely discussed || : 

3 

<^73> = -WiLrfp) = A E ^ K e ijK ■ (1.5) 

K=l 

The condensate Ql.5| ) breaks the original symmetry group SU(3) C (8) SU(3)l <S> 
SU(3) r (g> U(l) B down to 

SU(3) C+L+R ®Z 2 . (1.6) 



1 The presence of an external chromagnetic field may affect CSC: for an analysis see J_8| 
2 The more general case is 

(ipialpfp) = -(ipiai'jd) = 71 SaiSbj + 72 S a j5 b i , 

but the solution with 71 = —72 corresponding to (ll^q) is energetically favored. 
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Both the chiral group and the color symmetry are broken but a diagonal SU (3) 
subgroup remains unbroken in a way that locks together color and flavor (Color- 
Flavor-Locking=CFL model). We have 17 broken generators; since there is a 
broken gauge group, 8 of these generators correspond to 8 longitudinal degrees 
of the gluons, because the gauge bosons acquire a mass; one has 9 NGBs or- 
ganized in an octet associated to the breaking of the flavor group and in a 
singlet associated to the breaking of the baryonic number. The effective theory 
describing the NGB for the CFL model has been studied in |2C]. 

The 2SC model might be applicable for intermediate chemical potential /i, 
while for very high fi the CFL phase should set in; therefore the QCD phase 
diagram should display different phases as schematically depicted in Fig. 1.1. 



Quark Gluon Plasma 




Chemical potential 
Fig. 1.1 Schematic view of the QCD phase diagram 



Similarly to the chiral symmetry breaking condensate ( |1.1|), the phase tran- 
sition and the non vanishing of the condensates (|l.3| ) or (|l.5| ) result from a 
mechanism analogous to the formation of an electron Cooper pair in a BCS su- 
perconductor. At T = the only QCD interactions are those involving fermions 
near the Fermi surface. Quarks inside the Fermi sphere cannot interact be- 
cause of the Pauli principle, unless the interactions involve large momentum 
exchanges. In this way the quarks can escape the Fermi surface, but these 
processes are disfavored, as large momentum transfers imply small couplings 
due to the asymptotic freedom property of QCD. Even though interactions of 
fermions near the Fermi surface involve momenta of the order of n, their ef- 
fects are not necessarily negligible. As a matter of fact, even a small attractive 
interaction between fermions near the Fermi surface and carrying opposite mo- 
menta can create an instability and give rise to coherent effects. This is what 
really happens Q, || and the result is the formation of a diquark condensate, 
as expressed by ( |1.3| ) or Ql.5| ). We stress again that the only relevant fermion 
degrees of freedom are therefore those near the Fermi surface. 

This property and the effective theory for the fields at the Fermi surface is 
discussed for electromagnetism and ordinary BCS superconductivity in 22 



23] where the Wilson Renormalization Group approach |24[ is used to prove 
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the Landau Theory of the Fermi liquid p5[ and to prove the marginal character 
of the BCS interaction. First works towards the implementation of these ideas 
in the context of color superconductivity can be found in [26, 27, 28 1. In [29] 
this approach was used to build up an effective theory for the CFL model 
based on the approximation of the neglect of the negative energy states. This 
results in a rather terse formalism resembling the heavy quark effective theory 
fjOf and displaying as a characteristic note the existence of a Fermi velocity 
superselection rule and effective velocity-dependent fermion fields. Due to these 
features, we refer to this effective lagrangian as the High Density Effective 
Theory (HDET). In §1 the 2SC model has been studied by the same formalism, 
while in |32|, |33|, |3j] this effective theory has been applied to a different color 
superconducting phase, the crystalline color superconducting phase [35, 36, 37, 



38, B9J, also called the LOFF phase from the authors of the original model of 



the crystalline electromagnetic superconductivity [ 40 1 . 

The aim of this paper is to review the developments in the description of 
Color Super Conductivity that are based on the effective lagrangian approach, 
in particular the HDET approach. While there are good reviews of the whole 
subject, the approach based on effective lagrangian has not been extensively 
discussed yet and this paper aims to do that. In the remaining part of this first 
section we will present an elementary introduction to some basic aspects of color 
superconductivity, most notably the Wilsonian version of the Landau theory of 
Fermi liquids and the derivation of the gap equation. The effective lagrangian 
approach for the two basic models (CFL and 2SC) is discussed extensively in 
the subsequent Section 2. In Section 3 we review some theoretical developments 
that in a way or in another are related to the main theme of this survey. In 
Section 4 we present a detailed discussion of the LOFF state of QCD. Finally 
in Section 5 we discuss astrophysical implications of color superconductivity in 
the context of the compact stars phenomenology. 



1.2 Color superconductivity 

To start with, we consider a non interacting fermion gas; its energy distribution 

K E > = e P{E-n) + 1 ' 

becomes, when the gas is degenerate (T — ► 0), the step function, see fig. 1.2: 

f(E)^9(n-E) . (1.7) 

The condition E < fi defines a domain in the momentum space whose boundary 
is the surface E(p) = ep = fj>, called Fermi surface; €f is called Fermi energy and 
the momentum pp satisfying E(pp) = fi is the Fermi momentum. Any quantum 
state inside the Fermi surface, i.e. E(p) < ep is occupied by a fermion, and 
the other ones are empty. For m = 0, E = p and the Fermi surface is spherical 
Px + Py + p\ = Pf = l^ 2 ■ This results applies also for non relativistic non 
interacting fermion (with pp = 2m[i). 



6 



In presence of interactions the situation is more complicated, because, in 
general / = f(f,p). At very large chemical potentials one can however take 
advantage of the asymptotic (large p) properties of QCD, which implies that at 
T = the only QCD interactions are those involving fermions near the Fermi 
surface, as we remarked already. We shall discuss in the sequel the role of the 
surviving interactions for the quarks at the Fermi surface that, even if small, 
can be not negligible. For the time being we observe that the chemical potential 
is related to the fermion density (number of fermions per unit of volume n by 



d 3 p 

{2-Khf 



where g is the degeneracy of the state. If f(r,~?p) does not depend on f this re- 
sults in a constant density n. Clearly the density increases with fx; in particular, 
for degenerate free fermions one gets 

^9 [ PF 2 , ^9V% n ^ 

n = WWJ pdp = ^- (L8) 

In a similar vein the energy density and the pressure are given by: 

e = */(W^' (L9) 

P = 9fT^ P -f (HO) 
respectively. For non interacting fermions the state we have described is the 



and 



Fig. 1.2 Fermi distribution f(E) at T = 0. f(E) = for E > p = e F , and 
f(E) = 1 for < E < ii = e F . 

ground-state, i.e. the vacuum. Let us now suppose that the fermions are quarks 
and let us discuss the origin of the color superconductivity; we show by a simple 
argument that the naive vacuum is unstable provided an attractive interaction 
exists. Let us consider a pair of non interacting fermions near the Fermi surface; 
their energy is ~ Let us substitute them with a bound state, i.e. a pair 
(a boson) with binding energy ei n t < 0. The fermion energy decreases by the 
amount —2//; on the other hand the boson energy increases by 2fi + ei n t, so that 
the total energy is decreased by the amount — | | . Therefore the formation of 
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a bound state is favored. This might appear a small thing, but it is not, as one 
can immediately realize since there is an infinite number of ways to do that. 
The Cooper pair has preferably total momentum and total spin zero, but there 
is an infinite number of possible pairs, corresponding to the different directions 
of the Fermi velocity, see fig. 1.3 

v F = ?L . (i.n) 




Fig. 1.3 Cooper pair of two fermions living at the Fermi surface with total 
momentum and total spin zero. 

For the phenomenon of the superconductivity to occur one therefore needs 
an attractive interaction between the fermions. In the case of electronic su- 
perconductivity in the metals the attraction is the effect of the interaction of 
the electrons with the phonon: the resulting interaction between the fermions 
is indeed attractive and provides the required binding energy of the electronic 
Cooper pair. In the case of quarks the situation is more favorable because 
the QCD interaction in the color antitriplet 3 (antisymmetric color channel) is 
attractive; therefore color superconductivity can occur even in absence of a lat- 
tice. The result of the Cooper interaction is the formation of a colored diquark 
condensate, with total baryonic number 2. 

The possibility of color superconductivity was studied long ago ^ , but only 
recently § it has been realized that the effect offers several physical possibil- 
ities and is indeed sufficiently robust to survive the thermal fluctuations. The 
physical mechanism at its basis arises from interactions between two fermions 
at the Fermi surface. They will be discussed in the next paragraph. 

1.3 Effective action a la Wilson 

The previous results can be expressed in a more formal way by making use of 



the Landau theory of the Fermi liquid [25]; we will describe it in the Renor- 
malization Group language a la Wilson |l|' H H H 5 The idea of this 
approach is that the effective interactions of the fermions at T = can be de- 
scribed substituting the particles with their excitations, the quasiparticles that 
are in one-to-one correspondence with the particles. If one is interested only in 



3 Our presentation follows a similar discussion is in 
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special values of the momenta (in our case momenta near the Fermi surface), 
one can ignore the details of the interaction, which considerably simplifies the 
description without loss of physical insight. The simplification consists in the 
fact that the quasiparticles can be treated as free. However the form of the 
energy dependence on the momentum, i.e. the dispersion law e = e(p), differs 
from the free particle case. 

To derive in the framework of the Landau theory the results already dis- 
cussed in previous sections, let us start with a gas of fermions at T ~ and the 
free lagrangian 

£ = ip(i$ + Wo)ip , (1-12) 
from which the free action can be constructed 

S = J dtdp^(p){id t + n-e(p))i>(p) . (1.13) 

For free particles one obtains ( |1.13j ) using the Dirac equation to write a ■ pip = 
707 • pip = e(p)ip (see below, Eq. ( |2.lD ); therefore, for free particles, e(p) = 



\Jp 2 + m 2 . However the interactions can modify this simple dispersion law and 
therefore we do not use this result. Because of this possible modifications, as 
we have already stressed, the particles are substituted, in the Landau theory 
of the normal Fermi liquid, by quasi-particles, that we can imagine as particles 
dressed by their interactions. 

Let us now consider what happens when we scale energies and momenta 
down towards the Fermi surface. Let us divide the fermion momentum as 
follows 

p = k + i, (1.14) 

with k on the Fermi surface and e(k) = [i (i.e. k is a Fermi momentum), while 
£ = in, with ft a unit vector orthogonal to the Fermi surface. We scale down 
the momenta by a factor s — > so that they approach the Fermi surface; clearly 

E -> sE , 

£ -> si. (1.15) 
We note that these conditions can be also expressed as follows: 

E, \£\<S . (1.16) 

with 5 < fi. Moreover 

de 

dp dp 

where 



e(p) - fi -» e(k) - [i + (p - k) — ^ = i ■ — ^ = iv F {k) , (1.17) 

p=k 



Op 

is a vector orthogonal to the Fermi surface. Therefore 



v F (k) = || ^ (1.18) 

p=k 



S= I dtdkdi^ip) [id t - iv F (k)^j ip(p) . (1.19) 
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Under the scale transformation (1.15) we also have 



dt -> s^dt , d 3 p = d 2 kd£ -> sd 2 kd£ , 

d t ^ sd t , £^s£ , (1.20) 

which means that the fields ip in the action S of Eq. ( |1.13| ) scale as s~ 1 / 2 , 
while S scales as s . Let us note explicitly that for spherical Fermi surfaces the 
integration is 2-dimensional and the measure is 

dk=p 2 F dn. (1.21) 

Let us now consider the terms in the action that might result from the 
interactions; they are constrained by the symmetries of the problem. We classify 
them according to the number of fermion fields: 

1. A term quadratic in the fermion fields, such as 

dtd 2 kd£m(k)^(p)^(p) , (1.22) 



scales as s . Potentially this is a relevant operator PL However m(p) 
can be absorbed f\ in the definition of e(p), so that fll.22j ) is not a new 
operator at all. As to other quadratic terms such as, for example, those 
containing time derivatives and/or factor of £ they can be either absorbed 
in the corresponding terms in (|1.13j ) or are irrelevant. 

2. A quartic term such as 

dt uU ( d % di j) (^HpMps)) (frntipi)] 

x F({%})<5(pi + p 2 - p 3 - p 4 ) (1.23) 

scales as s = s _1 s _4//2+4 multiplied by the scale factor of the Dirac delta; 
it is therefore irrelevant if one can assume that 

S(pi +P2 - P3 - Pi) = S(k 1 +£i + k 2 + £2-k 3 -h-ki-£i) 

-» 5{h +k 2 -h~ ki) ■ (1-24) 



3. Operators containing 2n (n > 2) fermion fields scale as s n 1 and are 
therefore irrelevant. 

4 In the language of the Renormalization Group, an operator that becomes less and less 
important when the momenta scale is called irrelevant; an operator that becomes more and 
more important is called relevant; an operator equally important at all scales is called marginal. 

5 Since we can absorb m(p) into e(p) we cannot assume for the dispersion law the free particle 
expression; the specific form of the dispersion law depends on the nature of the interaction. 
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This is the Landau theory of the normal Fermi liquid derived by the renormal- 
ization group approach. Let us now explain how this theory is modified for 
superconducting materials. Let us consider an interaction term with a quartic 
coupling ( 1.23j ) such that the two incoming momenta p\ and P2 satisfy P1+P2 = 
(see Fig. 1.3). If the Fermi surface is parity symmetric, which we assume, then 
k\ + k2 = and also k^+k^ = 0. In other terms, for the special case of scattering 
with zero momentum, instead of flL24Q we have 



5(pi +P2-P3 ~Pa) 



(1.25) 



which scales as s . 

It follows that the corresponding operator becomes marginal and, since there 
are no relevant interactions, it may dominate. This is the mechanism at the 
basis of the phenomenon of Cooper pairing in the BCS f\ superconductors. The 
Cooper pairing dominates over the electronic electrostatic repulsion, in spite 
of the fact that the Cooper pair has a size of 10 4 A and the binding energy 
between the electrons is of order 10 -3 eV, while typical electronic interactions 
are of the order of the electron volt. 



-p' , E + E' 





p' , E-E' 

Fig. 1.4 Four-fermion vertex with the marginal coupling (sum of incoming mo- 
menta p+ {—p) = 0), together with the loop correction; +q and — q are the final 
momenta. 



The origin of the attractive interaction lies in the electron-phonon interac- 
tions in the metal. In order to study them we substitute them with an effective 
four-fermion interaction that, as we know, for special values of the momenta is 
marginal. 

The importance of the marginal interactions depends on the quantum cor- 
rections because they can produce the effect of reducing or increasing the cou- 
pling. Let us therefore show that, if one considers the loop diagram of Fig 1.4, 
the attractive interaction becomes stronger. 

6 The BCS theory, valid in the limit of weak electron-phonon coupling, is generalized to 
strong coupling by the Eliashberg theory MlJ. 
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To compute Fig 1.4 we use the coupling (1.23) assuming, for simplicity a 
constant interacting potential V. The fermion propagator can be easily derived 
from the ( |i~l~9| ). The energy and the longitudinal momentum integration are 
cut-off, so that one is actually considering only momenta near the Fermi surface. 

The sum of the two diagrams is given by 



V(E) = V - V 2 x 



I 



dE'd 2 kd£ 



(2tt) 4 [(£ + E')(l + ie) - v F (k)£][(E - E')(l + ie) - v F {k)l] 
5 
E 



V-V 2 Nln-, (1.26) 



N= I (1.27) 

(2TT) 3 V F (k) 



where 



and 5 is the upper cutoff in the energy integration. From ( |1.2q ) we get, at the 
same order of approximation, 

V(E) « T (1.28) 

which shows that, when E — > 0, if V > (repulsive interaction), V(E) becomes 
weaker and weaker while, if V < (attractive interaction), the coupling becomes 
more negative and stronger. Therefore one expects that, if V < 0, at zero 
temperature the naive vacuum is unstable and the formation of fermion pairs 
is favored. 

For metals at low temperatures the situation is as follows. If at some in- 

1 /2 

termediate scale E\ = (-^j) 5 (m = electron mass, M = nucleus mass) the 
attractive interaction is stronger than the Coulomb interaction, then at smaller 
energies the ordinary vacuum, filled with gapless fermions, will become unsta- 
ble and the electrons will condense; if, one other hand, this condition is not 
satisfied, the metal is an ordinary conductor. In a superconductor there is a 
non-vanishing vacuum expectation value of the electron bilinear: 

ftMfe) / . (1.29) 



The condensate ( |1,29 ) acts as an order parameter and its non- vanishing value 



signals the phenomenon of spontaneous symmetry breaking. As the broken sym- 
metry is a gauge symmetry, £7(l) e . m ., there are no Goldstone bosons associated 
to it, but the photon acquires a mass, which produces, for example, the Meiss- 
ner effect. This is what happens for BCS electromagnetic superconductors. For 
QCD with light quarks at zero density and zero temperature this mechanism is 
supposed to favor the formation of a qq condensate, with breaking of the chiral 
symmetry and associated Nambu-Goldstone bosons (the pions). For QCD at 
high density an attractive coupling (V < 0) between a quark pair is provided by 
the color antitriplet channel; therefore it should ultimately dominate for T — > 
and ji — > co. 
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1.4 The true vacuum 



Once we have been convinced that the naive vacuum is unstable we should 
try to find the true superconducting vacuum state. To do that we have to 
show that the formation of a quark-quark condensate occurs because of the 
color interaction. The interaction term is fundamentally provided by gluon 
exchanges and to compute it and derive the gap equation one can follow different 
approaches, such as, for example, the one of Bailin and Love in [|| based on the 
resummation of bubble diagrams. For this method and similar approaches we 
refer to the reviews J14|| . 

For pedagogical purposes it is however sufficient to consider a Nambu-Jona 
Lasinio (4^] (NJL) interaction (4 fermion interaction) that mimics the QCD 
interaction. In this approximation the hamiltonian is: 

H = J d 3 x (4>(i $ + /i7o)V> + F^T A ^j^T A ij) . (1.30) 

Here F is a form factor that is introduced to take into account the asymptotic 
freedom property of QCD; its Fourier transform behaves as follows: F{p) — > 
for p — > oo, and this property also ensures the good behavior of the integrals 
for p — > oo. 

In order to obtain the gap equation we make a mean field ansatz: 

i&C^pj - (^CVto). (1-31) 

We repeat that Greek indices a, (5 are color indices, Latin indices i, j are flavor 
indices. Let us now show that one expects 

(^aiC^fSj) = ^api^ijK^IK ■ (1-32) 



It is useful to note that, by evaluating ( 1.30| ) in mean field approximation, i.e. 



substituting to a pair of fermion fields in ( |1.30| ) their mean value given by ( 1.32j ) 
and ( p.. 31] ), the interaction term takes the form 

~ Aipip + h.c. , (1.33) 

i.e. the fermions acquire a Majorana mass. 



To prove ( 1.32 ) we note that, as we have already discussed, QCD interaction 
favors antisymmetry in color; the pair has zero angular momentum (in this way 
the entire Fermi surface is available and the effect becomes macroscopic); if it 
has also spin zero the fermion pair must be in an antisymmetric state of flavor 
to produce an antisymmetric wavefunction Q 

Clearly the form of the matrix £Iik corresponding to the true vacuum state 
depends on dynamical effects. In the literature two cases have been widely 
discussed 



7 With one flavor this possibility does not exist and one remains with the sole possibility of 
a condensate with angular momentum equal to 1; in this case, however the superconducting 
effect is more modest, with gaps of the order of 1 MeV |43|. 
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1. Color-Flavor-Locking (CFL): Q/x = 5ik, see ( |1 . 5| ) ; 

2. 2 flavor SuperConductivity (2SC): Qik = 8136x3, see (|Q|), 

and we will concentrate our attention on them (other cases have been also 
considered; we refer to the literature, as reviewed for example in fl4| ) 
The CFL model takes its name from the well known property 

3 

£a/3K£ijK = 8 a i5/3j — 5 a j5/3i , (1-34) 

K=l 

by which the color and flavor indices become mixed and indeed locked together. 

The 2SC model corresponds to the presence in the condensate of only the u 
and d flavors: this means that this ansatz may be relevant in the situation char- 
acterized by an intermediate chemical potential m u , <C fj, <C m s . Therefore 
one can imagine of having taken the limit m s — > oo and the theory is effectively 
a two-flavor theory. 

Let us now introduce a, a', b, ... annihilation/creation operators of particles 
and holes: 

i>ai = V u k (a k e lkx + b{e- lkx ) . (1.35) 

k 

The existence of the condensate gives rise to the following hamiltonian 
H = yi(\k- ulajak + (k + jj)b\b^j 

k 

+ ^^e^a k a- k + ^fe^blbl k + hc 

k 

where we have omitted for simplicity any color-flavor indices as we wish to stress 
the mechanism that produces the condensate. Let us now perform a unitary 
Bogolubov-Valatin pi| | transformation. To do that we introduce annihilation 
and creation operators for quasiparticle and quasiholes 

y k = cos 9a k - e l ® sin 6a^ k , 

cos ipb k — e** sin ipb^_ , . (1.36) 



By an appropriate choice of the parameters of the transformation we can trans- 
form the original hamiltonian in a new hamiltonian describing a gas of non 
interacting quasiparticles 

H = Y1 [ e y( k )ylyk + e*(fc)4**] (i-37) 

with 

€ y {k) = y/(k-li) 2 +A2F* , 



e z (k) = y/(k - m) 2 + A 2 F 4 . (1.38) 
These equations show two effects: 
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1. The quasiparticles and quasiholes are free, there is no interaction term; 



2. in the dispersion laws of the quasiparticles a mass term proportional to 
A, called gap parameter, appears. 

To obtain (|L37| ), i.e. a free hamiltonian, the parameters have to be chosen as 
follows: 

cos 26 



y/\k — /u| 2 + A 2 F 4 ' 

k + u , 
cos2oj = p (1.39) 

VI(fc + /i) 2 + A 2 F 4 

We note that, differently from the original operators that annihilate the 

false vacuum, the quasiparticles annihilation operators destroy the true vacuum 

y fc |0>=0, z fc |0>=0. (1.40) 



We have still to prove that A / 0. We do that and get an equation for A 
by substituting for y^, in eq. ( 1.32j ). We get in this way an integral equation 
that has the following schematic form (gap equation): 



A = C / dkF 2 (k) 



- A A 



;i.4i) 



We see immediately that the origin of the instability of the false vacuum lies 
in the first of the two terms in the r.h.s. of ( 1.41) ) : If A = there is no 
compensation for the divergence at k = fM. 

To be more quantitative we have to consider the two models in more detail. 
We differ this quantitative analysis to Section while giving here only the 
results. 



CFL model (m s = m u = = 0). 



In this case all the 3x3 = 9 quarks acquire a Majorana mass. The CFL 
condition gives two different set of eigenvalues. The first one comprises 8 degen- 
erate masses Ai = A2 = ... = As and the second set the non degenerate mass 
A9. The actual values of the gaps depends on the model and the approximations 
involved; typical values are 

Ai = ... = A 8 = A w 80MeV , 



A 9 « -2A , (1.42) 

for \i ~ 400 MeV (a numerical analysis for this case is contained in Section 
|2.6| ). The interaction couples quarks of different flavor and color. Moreover, 
since we start with a massless theory, the pairing occurs for left-handed and 
right handed quarks separately. Therefore these condensates break SU (3) c <8> 
SU(3)l<S>SU(3)r; finally, since the pair has baryonic number 2/3, also U(1)b is 
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broken. It should be observed, however, that a diagonal 5C/(3) c +l+/j symmetry 
remains unbroken. 



2SC model (m s >/j> m u ,m,i). 

In this case there are 4 massive quarks (the up and down quarks with colors 
1 and 2), while the strange quark with any color and the u and d quarks with 
color 3 remain ungapped. Numerical results for this case are qualitatively in 
agreement with the results ( 1.42j ) of the CFL model (see for example Section 



2J|. 

One can show, using the Schwinger-Dyson equation (also called in this con- 
text the Nambu-Gorkov equation), that, for arbitrary large (M, the CFL model 
is favored ||; however, for intermediate /x the 2SC state can be more stable and 
the phase diagram should resemble that depicted in Fig. 1.1. 

The spontaneous breaking of global symmetries implies the existence of 
Nambu Goldstone bosons; for internal symmetries there are as many NGBs as 
there are broken generators G a . The NGBs are massless scalar particles that 
interpolate among different degenerate vacua 

[H, G a ] = G a \0 > + . (1.43) 

As they are massless, they are the lowest energy quasiparticles of the effective 
theory. On the other hand, for broken local gauge symmetries, by the Higgs- 
Anderson mechanism the gauge bosons acquire masses and there are no NGBs. 
In CFL model the diquark condensate ([□]) breaks 577(3) L x SU(3) R x SU(3) C x 
U(1)b to SU(3) c+ l+r x Z2- All the 9 quarks are massive and they belong to 
a SU(3) c+ l+r singlet and a SU(3) c+ l+r octet. All the 8 gluons are massive 
and are degenerate. There are 8+1 Nambu-Goldstone bosons. In the 2SC 
model the condensate flT|) breaks SU(3) C SU(2) L ® SU(2) R x U(l) B down 
to SU(2) C x SU(2) L x SU(2) R . While the chiral group is unbroken, SU{3) C is 
broken to SU(2) C . Therefore 3 gluons remain massless and 5 acquire a mass. 
As to the other quasiparticles, there is one would-be NGB associated to the 
breaking of the axial color; moreover of the 6 quarks (2 flavors in 3 colors) 4 
are massive and 2 are massless. 



1.5 Effective lagrangian for the NGB in the CFL 
model 

We wish now derive an effective lagrangian for the low energy excitations 
(NGBs) of the CFL model p0|] . In order to do that let us introduce two coset 
fields X and Y transforming under the symmetry group as the left handed and 
right handed quark fields 

X - g c XUl , 

Y -> g c YUl. (1.44) 
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X and Y contain both 8 scalar fields 

X = exp i.i 

Y = exp < i 



f[ A T A 
F 

fl A T A 
F 



(1.45) 



corresponding to the 16 NGBs induced by the spontaneous symmetry breaking 
SU(3) C x SU(3) L x SU(3) R -> SU(3) C+L+R (we neglect for the moment the 
U(1)b factor). We introduce the gauge covariant currents 



J 1 ' 

Jy 



Yd»Y ] + f , 



(1.46) 



where 

9fl = igT A A A (1.47) 
is the gluon field. Under the gauge group these currents transform as follows: 



J x 



Jy 



3cJ x 9, 



9cJy9t ■ 



(1.48) 



The most general gauge invariant (but not Lorentz invariant) lagrangian con- 
taining X, Y and g^ fields is as follows 



Ceff = -^9,u{Tr[(J^-4)(r x - 

+ aTr[(J x + j£)(J x + JZ)}} 
+ (kinetic term for g) , 



y) 



(1.49) 



where g^ u = diag(l, —v, —v, —v) and F, a are parameters. We will prove 
in par. 2.3.2| that an effective lagrangian of this type for the scalar fields is 
generated by loop corrections. We rewrite this equation as follows 



C 



eff 



4 

+ aTr 



g^jTr \{Xd^ - Yd^){Xd y X^ - Yd y Y^) 



{X8^ + Yc^yt + 2 9tl )(Xd v Xi + Y3 V Y^ + 2<7„)] } 



+ (kinetic term for g) . 



(1.50) 



Note that in this equation one has both a Debye and a Meissner mass for the 
gluon, m? D = ag 2 F 2 , m 2 M = v 2 am 2 ~ ) corresponding to a mass term for Aq and 



A, respectively; they can be read from ( p.. 50 ). We will show in par. 2.3.2 that 
one has to put a = 1 in order to have the correct normalization for the gluon 
masses. Since F = O(n), in the limit of small momenta, one can neglect the 
kinetic gluon term and, because of the equations of motion, the gluon field 
becomes an auxiliary field given by 



9fi 



1 



{Xdprt - YduY^) ■ 



(1.51) 
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Therefore the lagrangian (1.50) reduces to 

C*ff = - ^9^Tr [{XdpXi - Yd^){Xd v X^ - Yd u Y^) 



If we put 

£ = Y ] X = exp 

T A = X A /2, from flL5p we finally obtain 

F 2 



2 iU A T A 
F 



eff 



Tr ( anas 1 " - v 2 vs • vs f 



(1.52) 
(1.53) 

(1.54) 
(1.55) 

The values of the parameters in (|1.52j ) have been derived by several authors 
(the derivation in the framework of the effective theory will be presented in 
Section |). One finds g|, HI H EH> U H : 

u 2 (21-81n2) 2 1 

' W = 3 • 



We note that the matrix £ transforms under the chiral group as 



F 



36vr 2 



(1.56) 



A similar analysis can be done for the superfluid mode [20],[45|; we can 
therefore write down the low energy effective lagrangian including both the 
SU(3) and the U(l) NGBs as follows: 

,? 2 

£-NGB 

+ kr (d t Ud t U ] - \VU ■ VU^ , (1.57) 



where 



U 



(1.58) 



sec below, eq. (|2~75| ). 

The lagrangian ( |1.49 ) is completely analogous to the hidden gauge symme- 
try version of the chiral lagrangian of QCD |51| |, but for the breaking of the 
Lorentz invariance: g^ u / g^ v , due to the high density. 

We finally note that there is an extra U(1)a symmetry in the lagrangian, 
which is broken by instantons. Even though the contribution of the instantons 
is parametrically small at asymptotic densities, it is sufficient to forbid the 
presence of an extra massless scalar 0. Notice that one has no NGB associated to 
the breaking of the chiral symmetry, because, as the non-vanishing condensates 
only couple quarks with the same helicities, chiral symmetry is broken only 
indirectly, by the fact that the mechanism of Color Flavor Locking locks left 
handed quarks to color as well as right handed quarks to color. 



Nevertheless, because of this suppression, this extra pseudoscalar is expected to be rather 
light. It can be observed that the instanton interaction is suppressed not only by the asymp- 
totic energies, but also by the fact that, being a 6-quark coupling, it becomes irrelevant at the 
Fermi surface. 
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Chapter 2 

High Density Effective Theory 



2.1 The 2-D effective theory 

We wish now to derive an effective description of QCD at very high baryonic 
chemical potential fi and very small temperature (T — ► 0). 

To begin with we consider the equations of motion for a massless quark at 
finite density 

(ii? + ^70) ip = , 

where = is the quark number chemical potential. The ^70 term corresponds 
to a /jlN = fii/j'i/j term in the hamiltonian. We note explicitly the breaking of 
Lorentz invariance of this equation, which, however, is still symmetric under 
rotations. At T = the vacuum is characterized by the fact that the fermions 
fill in all the low energy states up to a maximum energy, the Fermi energy. Due 
to the Pauli principle, the interaction of low energy quarks implies the exchange 
of high momenta, which is disfavored by the Asymptotic Freedom property of 
QCD. Therefore only the fermions that are near the Fermi surface can interact. 
Since they have high momentum (p ~ fi) they are almost free and we can write: 

i.e. in momentum space: 

ii> + ^70) i> = , 

or (a = 707): 

p ^ = (a -p- fi)ip . (2.1) 
Let us now consider the eigenvalue equation for the Dirac hamiltonian: 

H D = po = a ■ p - fi , (2.2) 

whose eigenvalues can be found solving the equation (po + m) 2 = P 2 - There are 
positive and negative energies, given by 

p = E ± = -fi± \p\ (2.3) 
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with eigenvectors 

V>± = P±^ = \ (l ± ^p) i> ■ (2-4) 

Therefore at high-density the energy spectrum of a massless fermion is described 
by states |±) with energies given by (^). As we have already stressed, for 
energies much lower than the Fermi energy /i, among the states |+) only those 
close to the Fermi surface, i.e. with \p\ ~ fi, can be excited. As to the 
states |— ), they have E_ pa —2p» and therefore decouple; this means that in the 
physical spectrum only the positive energy states |+) and the gluons remain. 

After this quantum mechanical introduction, let us consider the field theo- 
retical version of the previous argument. The main idea of the effective theory is 
the observation that the quarks participating in the dynamics have large (~ //) 
momenta. Wherefore one can introduce velocity dependent fields by extracting 
the large part \xv of this momentum. One starts with the Fourier decomposition 
of the quark field tp(x): 

and introduces the quark velocity by 

= ^ + p , (2.6) 



where 

with \v\ = 1. Let us put 

and 

with 



= (0, v) 

i= vi\\ + i± 



£ ± =£-(£■ v)v . 
We can always choose the velocity parallel to p, so that £± = and 



(fp = n* I <m I d£\\ I d£ = 4vr^ z / — / d£\\ \ d£ . (2.7) 



DC 



In general we will be interested only in the degrees of freedom at the Fermi 
surfaceR therefore we can limit the integration limits as follows: 

S < £\\ < +5 

with 5 < fi and 5 much larger than any other energy scale, in particular the 
energy gap masses, if any. By this hypothesis 

p+S 

d 3 p = 4vr / u 2 ^ ' / d£\\ , (2.8) 



1 There are some important cases where this limitation does not apply, see below. 
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with 

V= f — 

V 

eventually we take the limit 5 — > +00. In this way the Fourier decomposition 
(p.5|) takes the form 



where 

d 4 £ 4tt^ 2 



w - w • dH ■ <2 ' 9) 

while are velocity-dependent fields. The projection operators in (2.4) can 

be approximated as follows: 

ft = ^ (2.10) 

and one can write 

Mx) =^2'e-^ W + (x)+il>-(x)] , (2.11) 

v 

where 

Mx) = P±Mx) =P±j ^e-^A<i) ■ (2-12) 

It is clear that ip± are velocity dependent fields corresponding to positive and 
energy solutions of the Dirac equation. 
Let us now define 

V» = (1, tO , V» = (1, -v) , 

7f = (7°, (v^)v), 7 ^= 7 ^- 7 f. (2.13) 

Moreover we use 

^_7^V+ = 1^-7^+ • (2-14) 

Substituting into the Dirac part of the QCD lagrangian we obtain, in the hy- 
pothesis that the residual momenta £ are small, 

+ (^+ii? ± ^_ +h.c.)l ; (2.15) 

= D^jj^ and D^, is the covariant derivative: -D M = + igA^ . We note that 
here quark fields are evaluated at the same Fermi velocity; off-diagonal terms 
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are subleading due to the Riemann-Lebesgue lemma, as they are cancelled by 
the rapid oscillations of the exponential factor in the fj, — > oo limit. One may 
call this phenomenon Fermi velocity superselection rule, in analogy with the 
behaviour of QCD in the uiq — > oo limit, where the corresponding effective 
theory, the Heavy Quark Effective Theory exhibits a similar phenomenon f30[| . 
By the same analogy we may refer to the present effective theory as High Density 
Effective Theory (HDET), 

We can get rid of the negative energy solutions by integrating out the -in- 
fields in the generating functional; at the tree level this corresponds to solve 
the equations of motion, which gives 

iV-Dip + = (2.16) 

and 

V>_ = % -^— 7o^+ • (2.17) 

2 (j, + iV • D 

This equation shows the decoupling of ^_ in the fi — > oo limit. Therefore, in 
this limit, the ip- field plays no role. We also note that in the equation for tp + 
only the energy and the momentum parallel to the Fermi velocity are relevant 
and the effective theory is two-dimensional. 

We can write the effective lagrangian at the next to leading order in fi as 
follows: 

1 



£ = --F« u F a ^ + £ D 



where 



4^ L 2fi + iV-D 



(2.18) 



Using the identity 
where 



y 2 L 



V^V" + V V V' M 



we can write (|2.18|) as follows 



if>Uv ■ Dtp - p^v 1 " 



1 



2fi + iV ■ D 



(2.19) 



(2.20) 



We can make this equation more symmetric by introducing positive energy 
fields with opposite velocity: 



(2.21) 



Even though we use the same notation ip± here and for the positive/negative 
energy fields, we stress that ip± in ( |2.21| ) and in the sequel are positive energy 
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fields, with opposite Fermi velocities: iv. Eq. ( p. 20 ) can be therefore written 
as follows: 



Co 



E K 



iV-Dtp + + ipUv-Dip 



pfiu 



'2/1 + iV ■ D 



Z? M A,^ + - ip 



t 



2ji + iV ■ D 



. (2.22) 



To this lagrangian involving only the left handed fields, one should add a similar 
one containing the right handed fermionic fields. We observe that a factor — 
has been introduced in the definition of the sum over the Fermi velocities that 
now reads: 



E 



dv 
8^ 



(2.23) 



This extra factor is needed because, after the introduction of the fields with 
opposite v, ip- = ip + _tf, we have doubled the degrees of freedom and therefore 
we must integrate only over half solid angle. Apart from this factor there 

should be another ^] factor - in ( 2.22| ), but we got rid of it by a redefinition of 

the fermion field, 







V2 



We close this section by a definition that will be useful in the sequel: 

C D =C + C 1 + C 2 , (2.24) 



with 



Co 



V 

i g Y yP+iV ■ Aip+ + ipUV ■ Aip. 



(2.25) 
(2.26) 



and 
£2 







1 



2/1 + iV D 



1 



2fj, + iV ■ D 



D„D u i;- . (2.27) 



Co is the free quark lagrangian (without gap mass terms); C\ represents the 
coupling of quarks to one gluon; £2 is a non local lagrangian containing cou- 
plings of two quarks to any number of gluons. Its effect will be discussed in 
section 



2.2 HDET for superconducting phases 

The construction described above is valid for any theory describing massless 
fermions at high density provided one excludes degrees of freedom very far from 

2 In |29|] this factor explicitly appears in the lagrangian; here we have simplified the notation. 
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the Fermi surface. As discussed in Section [j], however, for small temperature 
and high density the fermions are likely to be gapped due to the phenomenon of 
the color superconductivity, We shall examine the modification in the formalism 
in two models, i.e. the Color Flavor Locking, CFL, model (3 massless quarks: 
m u , m<2, m s <C /i) and the Two-Flavor, 2SC, model (m s — ► oo) 



2.2.1 HDET in the CFL phase 

In the Color Flavor Locking phase the symmetry breaking is induced by the 
condensates 

($ L JC^%) = -^C^%) = | e^jeyj , (2.28) 

where ip ,R are Weyl fermions and C = i&2- Eq. ( p. 28 ) corresponds to the 
invariant coupling (ip = ipi)'- 

"f 2 VCei1>ei -{L^R) + h.c, 

7=1,3 

and 

(ei) ab = e Iab . (2.29) 
For the Dirac fermions ip± we introduce the compact notation 

*-(<&) {2X) 

and we use a different basis for quark fields: 

1 9 

V>± = -7= E A ^± • ( 2 - 31 ) 

The CFL fermionic lagrangian has therefore the form: 

C-D = £q + £-1 + £a = 

E V a\ ( iTr[T A V ■ DT B ] A AB \ B 

x { A AB iTr[T A V-D*T B ]) X 

+ (L -> R) (2.32) 

where 

A AB = ATr[ejT A ejT B ] (2.33) 

and 

Ta = ^- (2.34) 
Here Ag = Ao = x 1. We use the identity (g any 3x3 matrix): 

e I g T e I = g - Trg ■ (2.35) 



24 



we obtain 

A AB = A a 8ab (2.36) 

where 

Ai = ■ ■ ■ = A 8 = A (2.37) 

and 

A 9 = -2A . (2.38) 
The CFL free fermionic lagrangian assumes therefore the form: 

E E i^o) x A + (^R)- (2.39) 

v A=l V 7 

Here the average over the velocities is given by (|2.23j ). From this equation one 
can immediately obtain the free fermion propagator that in momentum space 
is given by: 

Sab(p) = - t ~/ A * , . ( V 7 ~ A t) (2-40) 



V ■ p V ■ p - A\ + ie \~ Aa V ■ p 

We note that this propagator is analogous to the Nambu-Gorkov propagator of 
the full theory. 

2.2.2 HDET in the 2SC phase 

For the two flavour case, which encompasses both the 2SC model and the ex- 
isting calculation in the LOFF phase we follow a similar approach. 
The symmetry breaking is induced by the condensates 

corresponding to the invariant coupling (ifj = ip L ): 

C A = -— ^ T Ceipe - (L -» R) + h.c. , 

where 

e = ia 2 • (2.41) 

As in the previous section we use a different basis for the fermion fields by 
writing 

^ + ,ai = (<,a = l,2) 

A=0 V 2 

V>+,32 = r + , (2.42) 

where a a are the Pauli matrices for A = 1,2,3 and uq = 1. As usual ^+ are 
positive energy, velocity dependent fields and we also introduce V>_ according 
to the formula: 

^ = <±* ■ (2-43) 
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A different, but also useful notation for the fields tp+ a i uses a combination 
of A matrices, as follows 

^ = E ^itf • (2-44) 
The Xa matrices are defined in terms of the usual A matrices as follows: Ao = 

^n^ 8 + V^ o; ^ A = ^ A ^ = 15 2 ' ^ 4 = y^ 5 ; ^ 5 = 

After the introduction, analogously to (|2,30| ), of the fields x" 4 ; the 2SC 
fermionic lagrangian assumes the form: 

£-D = A) + A + £a = 

V V v^t ( lTr ^ A V ' Df B] ^AB \ n 

^ ^ X \ A AB iTr[T A V ■ D*T b \I 



V 



A,B=0 



Here 



and 



+ (L -»■ R) . (2.45) 



A AB = ^Tr[ea T A ea B ] (A,B = 0,...3) 

A AB = (A 5 = 4, 5). (2.46) 



T A = ^| (A = 0,...,5). (2.47) 



Analogously to ( 2.35| ) we use the identity (5 any 2x2 matrix): 

eg T e = g-Trg; (2.48) 

we obtain 

A AB = A A 5 AB (2.49) 

where is defined as follows: 

A A = (- A, +A, +A, +A, 0, 0) . (2.50) 

Therefore the effective lagrangian for free quarks in the 2SC model can be 
written as follows 

From this equation one can immediately obtain the free fermion propagator 
that in momentum space is still given by ( [2.40 ), with the A^ given by ( 2.50) ). 
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2.3 NGB and their parameters 



In the CFL model the diquark condensate breaks SU (3) l <8) SU (3)r <8) SU (3) c <g> 
U(1)b to SU(3) c+ l + ji (g) Z 2 . By the Goldstone theorem we have 8 Nambu- 
Goldstone Bosons +1 related to the breaking of U(1)b- We have described the 



effective theory for this case in par. 1.5. On the other hand in the 2SC model 



there is only one would-be NGB associated to the axial color because the chiral 
group is unbroken and there is an unbroken subgroup that plays the role of 
U{l) B . 

How to describe the NGB couplings to the fermions? In order to derive an 
effective low energy lagrangian for the NGBs and obtain these couplings we use 
the gradient expansion (for a clear exposition see e.g. [jS^] ). First we introduce 
as many external fields as the NGBs, with the same quantum numbers of the 
goldstones; subsequently we perform a derivative expansion of the generating 
functional. This gives rise to the effective action for the NGBs. 

2.3.1 Gradient expansion for the U(l) NGB in the CFL model 
and in the 2SC model 

To describe the procedure of the gradient expansion we consider the effective 
lagrangian for the NGB associated to the spontaneous breaking of the U{\) 
symmetry. For pedagogical purposes this calculation will be described in detail. 

To start with we consider the CFL model. The method of gradient or 
derivative expansion proceeds as follows. Let us introduce the external field a 
by the substitution 



C^lJ \UC^*_ 
where 

U = e ialh . (2.53) 
The quark lagrangian ( |2.3S| ) becomes 

a, + c A + c. = E E x Ai (l v J> * A + < L - R > ■ < 2 - 54 > 

v A=l V 7 

At the lowest order in the field a we have 



v A=l 



( o — - — ^ 

fa fa 

2ia 2a 2 



X A + (L - R) , (2.55) 



V fa P. J 

which contains the couplings crxx an d craxx- Since a is an external field there 
is no associated kinetic term. However such term is generated by quantum 
corrections. To show this we consider the generating functional depending on 
the external fermion sources 77 and (we consider only left-handed fields for 
simplicity): 

W[T],rf} = j[dx][dx f ][da]expli^2jx<Ax + rfx + x'ri\ (2.56) 
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where we have introduced 



and 



J a Jo 



r °" 1 °i + o ) ' ri " (l 



(2.57) 
(2.58) 



Clearly we have 

W[ V , 7?t] = I [d x ][d X j ][da}x 



x expji^y (x ] + ri ] A~ 1 ^ A( X + A~ l rj) - A' 1 ^ . (2.59) 
After a change of variables in the functional integration we obtain 
W[ri, if] = J[dx}[dx^[da}exp^i^2 J X f 4x - j 
= J[da]expl-iJ2j ^A^vj f[dx][dx ] ] expLj^J X^x^j 
= J [da] exp ^— i ^ 77^^4~ 1 ?7 j- exp {Tr In ^4} 



[da] exp < i 



iTr ln^4 — t^t4 ^ 



(2.60) 



Now we have 



„_i ^ (-1) 2 " / 2iaA f-2W 2 A 
: — i Tr In 5 +i >_ ( iS — - — iFo + iS ^ iTi 



n=l 



n 



fa fa 

This is a loop expansion. At the lowest order it produces the effective action 

S eff 



+— Tr I dxdy 



iS(y,x)2ia(x)A i S(x,y)2ia(y)A 

— i r -, i r 



fo 



fo 



+ iTr / dx 



iS(x,x) (-2) Aa 2 (x) 
ft 



(2.61) 



The two terms correspond to the diagrams in Fig. 2.1, i.e. the self-energy, Fig. 
2.1a, and the tadpole, Fig. 2.1b. They can be computed by a set of Feynman 
rules to provide 

iCeff ■ (2.62) 
In momentum space the Feynman rules are as follows: 
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Fig. 2.1 One-loop diagrams. External lines represent the currents J®, J^. Full 
lines are fermion propagators. 



1. For each fermionic internal line with momentum p the propagator 

i &AB fV-p -A A 



iSab(p) = i§ABS(p) 



V -pV-p- A\ + ie \-A A V-p 



; (2.63) 



2. For each vertex a term i£i n t that can be derived from the effective la- 
grangian; for example the a couplings to quarks can be derived from 



S); 

3. For each internal momentum not constrained by the momentum conser- 
vation the factor 



(2tt)< 



471-3 



+6 



(2.64) 



4. A factor of 2 for each fermion loop to take into account the spin (L + R); 

5. A statistical factor arising from the Wick theorem if needed. 

The result of the calculation of the effective lagrangian in momentum space is 
as follows: 



v, A,B ' 



x Tr 



ib A B {£ + p) — -, — iToibBA {£) — -, — «r 



47T 3 



fa 

d 2 lTr 



v,A,B 



is AB (ey 



fa 

-2A B a 2 



■»ri 



(2.65) 



corresponding to the two diagrams of Figs. 2.1a and 2.1b respectively. After 
some computation one has 



'' c eff{p) = i£i(jp) + iCiiip) = ^2 



v. A 



x / d 2 £ 



V ■ (I + p)aV ■ I a + V ■ {I + p)aV ■ la - 2A\a 2 
D A (£ + p)D A (£) 
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2a 2 - 
where we have put 

Da(p) = V -pV -p- A\ + ie . 
One can immediately notice that 

£/(p = 0) + £n(p = 0) = . 



(2.66) 
(2.67) 
(2.68) 



This result implies that the scalar a particle has no mass, in agreement with 
Goldstone's theorem. To get the effective lagrangian in the CFL model at the 
lowest order in the a momentum we expand the function in p (\p\ -C | A) to get, 
in momentum space 



iC e ff{p) = 
where we have defined 



v, A 



(V ■ p)a(V ■ p)a I 2 , 



(2.69) 



(2.70) 



(for other integrals of this form see the appendix at the end of this section). 
Since we have 

—in 

h 



2A\ ' 



one gets, in the configuration space, 

9/x 2 



C e ff(x) 



^T(V-da)(V-da) . 



(2.71) 



(2.72) 



Since 



Yv^v" = - [ —v»v v = - 

V 2 J 4vr 2 



(I \ 


0-: 

\0 



i o o 



V 



we obtain 



with 



Ceff(x) 



9/i 2 1 
vr 2 / 2 2 



(<V) 2 - ^(V.) 
1 

" 3 ' 



(2.73) 
(2.74) 



This kinetic lagrangian has the canonical normalization factor provided 

(CFL) , (2.75) 



_ V 

J a o 



IT 



which fixes the a constant. Therefore the effective lagrangian for the NGB a 
particle is: 



C 



1 



(d a) 2 - l(Va) 



(2.76) 
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We note that the value of the velocity (2.74) is a consequence of the average 
over the Fermi velocities and reflects the number of the space dimensions, i.e. 
3. Therefore it is universal and we expect the same value in all the calculations 
of this type. 

In the case of the 2SC model the formalism applies with obvious changes. 
In particular again we introduce the external field a as in (2.53); the effective 
quark-NGB lagrangian is given by the analogue of ( |2.39 ) where now the sum 
over A runs from A = to A = 5 and is given by ( |2.5C|), T he final result 
differs from the CFL case only in the coefficient in front of ( p.75[) which reflects 
the number of color-flavor gapped degrees of freedom, 9 in the CFL case and 4 
in the 2SC case; therefore one has 



,2 = V 
7T Z 



(2SC) 



(2.77) 



whereas the result (|2.74f) being universal holds also in this case. The NGB 
effective lagrangian is still given by (2.76). The NGB boson is in this case 
only a would-be NGB because the axial color is explicitly broken, though this 
breaking is expected to be small at high densities. 



2.3.2 Gradient expansion for the 577(3) NGBs 

The gradient expansion for the SU(3) NGBs is more complicated because of 
the color-flavor structure of these fields. However, the procedure is conceptually 
identical to the case of U(l): the NGBs are introduced as external fields and 
acquire a kinetic term, thus becoming dynamical fields, by integrating out the 
fermion fields. Let us outline this calculation in the present section. 
From the invariant coupling 

-yE ^Ceiipe! (2.78) 

7=1,3 

the coupling of the left-handed Weyl spinors ip's to the octet of external scalar 
fields is: 

7=1,3 

with an analogous relation (with matrix Y) for the right handed fields. Both 
X and Y have v.e.v. given by 

(X) = (Y) = 1 

and shall use the gauge 

X = y f . (2.79) 

The lagrangian giving the the coupling of the quarks to the external fields can 
be obtained from (|2.32| ) and is given by 

£ D = EE x ^( iV : d5AB *ba ) x b 

v A,B=1 V Aa Aa/ 
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+ (L - R) , 



(2.80) 



where 



= A Tr 



T 



= A ^Tr[T A X^T B X^] - Tr[T A X ] \Tr[T B X^ . (2.81) 

One can now expand the NGB fields 

'A a IF 



X = exp i 



2F 



a = 1 



; 3 w : 



and obtain the 3-point xxH, and the 4-point xxHII couplings. The result is, 
for the 3-point coupling 



E z ^{E^[^r.x" + x»*r 0X » 

a=l V 



^xxn 

- £ ^X af r X 6 n c } . (2.82) 

a,6,c=l 

On the the other hand the 4-point coupling (2 Goldstones + 2 fermions) 
^xxnn 



a=l 



E 



+ 



a,fe,c=l 



£ ^2^x ct rox d n a n fe } 



(^2^x ct rix 9 n a n 6 + h. c .) 



(2.83) 



a,b,c,d=l 



~ 8 4 

Rafted = 2 £ (9cap9dbp + d c dpd a bp) ~ -Sac^db + -^cd^afc • 



p=l 



Using the Feynman rules given above and the interaction lagrangians ( [2.82 ) 

(2.84) 



and (2.85 we get the effective lagrangian as follows: 

km ^(21-81n2) l * / a a 1. \ 
C '"~ 36^ 2^V n 11 - 3 |VIIa| ) ■ 

a=l v J 

The canonical normalization for the kinetic term is obtained provided 

// 2 (21 - 8 In 2) 



F 



36vH 



(2.85) 



The NGB velocity is again l/\/3 because of the universality noted above. 
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2.4 Gluon masses and dispersion laws 



2.4.1 Evaluating the bare gluon mass 



Let us go back to Eq. ( p. 27 ). £2 is a non local lagrangian containing an infinite 



number of local terms; let us discuss in the formalism of the two component 
fermion field x introduced in eq. ( 2.30Q : 



X 



Let us work in the quark basis x = Xaii ( a =color, j =flavor); £2 can be 
written as follows: 

Let us consider this coupling in the momentum space where 

i-itu + gA c u \ c Il)(-i£ v + gA d v \ d /2) 

{2fi + iV ■ D\- 1 P tw D ll D v = — - t 1 A y v 1 ' P^ . (2.87) 

1 ; M 2^ + iF- (-i£ + ff A e A e /2) V ; 

Since the effective theory is bidimensional we can put £j_ = £ p P pfl = 0. The 
simplest coupling generated by this interaction term is a quadrilinear coupling 
(two quarks, two gluons); it is obtained by putting to zero the coupling to the 
gluon in the denominator and one gets: 

txxAA = -9 2 Y, X^TT i^f 1 -L) P " UX • (2 - 88) 
z z \ 2 fJ ,+v-eJ 

Let us now consider the CFL model. Introducing the quark basis x A (A = 
1, 9), we get the coupling in the following form: 

txxAA = ~9 2 £ x\A^At ( _jM P^Xfl^ . (2-89) 
where the Clebsch Gordan coefficient is given by 

£d* = ^tr(X A X c \ d X D ) . (2.90) 

£ XX /L4 contributes to the gluon self-energy through a tadpole graph analogous 
to that in Fig. 2.1b (the external lines should be interpreted now as gluon 



lines). It is important to stress that, due to the denominator in (2.8S), the 
loop integral takes relevant contribution from the region where £q ± £11 ~ — 2fi. 
This means that one cannot integrate only near the Fermi surface, but the 
longitudinal integration must be extended up to 2/i, or, in terms of £u = p — /x, 
from —fj, to +/U. 
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The tadpole contribution to the gluon self energy can be computed by using 
the set of Feynman rules discussed above and by ( 2.8S| ). One gets 



/+M (jPu r+oo 



d£ 
— x 

2-7T 



where D{£) = l\ - t\ - A 2 + ie and D'(jt) = l\ - £ 2 - A' 2 + ie. The factor 
4=2 x 2 originates from a statistical factor and a spin factor. Since 

E^ = k d ( 2 - 91 ) 



3 



and 



we obtain 



£ P"" = --S^S^ (2.92) 
3 



= o , = o , (2.93) 

D% = -wj^'tfcd (CFL) (2.94) 



where w p = g[i/\/2ir is the plasma frequency and we have used the result 
/ - ^2 f + °° d£ V-l in 2 , 



/' 



2vr v ^ 11 7_oo 2vr + V ■ £ + ie){£ 2 - £ 2 - A 2 + ie) 8n 



For the 2SC case the calculation is similar; one obtains the following results 
(c,d= 1,--- ,8): 

n°° = o , = o , (2.96) 

Kd = -™ 2 /i6 cd (2SC) (2.97) 

where m g = g/ji/ySTT. Note that the tadpole graphs only contributes to the 
Meissner mass. 

2.4.2 Debye and Meissner mass and gluon approximate disper- 
sion laws: CFL case 



For both the CFL and the 2SC model the calculation of Section |2.4.1 provides 



a mass term to the gluons; it is a chromo-magnetic mass term (Meissner mass) 

3 We note that in this case the tadpole diagram presents infrared divergences (in to). To 
control these divergences one considers the system in a heat bath and substitutes the energy 
euclidean integration I4 = — i£o with a sum over the Matsubara frequencies I4 — » uj n — 

2n ( n + — ] T; eventually one performs the limit T — * 0. For details see f3l| and the appendix. 
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since it gives mass to A; there is however a second diagram contributing to the 
Meissner mass, i.e. the graph of Fig. 2.1a, with external line interpreted as a 
gluon. As we shall see in this section it contributes also to the chromo-electric 
(Debye) mass and to the dispersion law, i.e. 

E = E(p) ; (2.98) 

masses and dispersion laws can be derived by the polarization tensor 

n^(p) . (2.99) 

For example 

n 00 (o) = m 2 D 

U ij (0) = -5 ij m 2 M . (2.100) 

To compute H^ v from Fig. 2.1a we need L\ i.e. the interaction term between 
the quarks and one gluon. The fermionic current could be extracted from eq. 
( p.26|). H owever it is more convenient to work in the quark basis \ A - We use 
eq. ( 2.32| ) to write, in the CFL case 

A = igA%J°, (2.101) 

where 



where 



We obtain 



where 



h AaB = -Tr[T A X a T B ] . (2.103) 



+ ^±^(v T _.oj„ ( , 104) 



9abc = d abc + if abc . (2.105) 



and d a bc, fabc are the usual SU (3) symbols. The result of the self energy diagram 
(Fig. 2.1a, with external lines representing fermionic currents) can be written 
as follows: 
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where the propagator is given by eq. ( p.63|) . We note the minus sign on the 
r.h.sof (12111 , due to the presence of a fermion loop and the factor 2 due to the 
spin (L + R). To this result one should add the contribution arising from the 
tadpole diagram of Fig. 2.1b (with, also in this case, external lines representing 



fermionic currents J", J*); this diagram was computed in 2.4.1: 



LL ab 



6 ab 5^5 M 5 vk iol . (2.107) 



lb p 



To derive the dispersion law for the gluons, we write the equations of motion 
for the gluon field A b ^ in momentum space and high-density limit: 

[(V/" + pV) <U + nj] 4=0. (2.108) 

We define the invariant quantities Ho,ni,n2 and H3 by means of the following 
equations, 

r n 00 = n (p , p) 
n^(p , p) = { n 0i = rr° = n^, P) n l (2.109) 

[ IPJ = U 2 (p , V) 5 ij + n 3 (p , V) n l ni 

with n = — . It is clear that in deriving the dispersion laws we cannot go be- 
p 

yond momenta at most of the order A, as the Fermi velocity superselection rule 
excludes gluon exchanges with very high momentum; it is therefore an approx- 
imation, but nevertheless a useful one, as in most cases hard gluon exchanges 
are strongly suppressed by the asymptotic freedom property of QCD. 
Besides the longitudinal and transverse gluon fields defined by 



A? L = (n-A a ) n t , 

A a T = A«-A<} L , (2.110) 
it is useful to consider also the scalar gluon field Aq. By the equation 

PuKtK = - (2- m ) 

one obtains the relation 

( Po n - \p\u 1 )A = n-A( Po n 1 - |pKn 2 + n 3 )) , (2.112) 

between the scalar and the longitudinal component of the gluon fields. The dis- 
persion laws for the scalar, longitudinal and transverse gluons are respectively 

(n 2 + n 3 + pi) (|p1 2 + n ) = Po \p\ (2^ + Po \p\) , 
(n 2 + n 3 +pI)(\p\po + n ) = Po |p| (2n x + P g) +nf , 

p 2 -\p\ 2 + U 2 = 0. (2.113) 

The analysis of these equations is rather complicated [^] . We give an approxi- 
mate evaluation of them by developing W IU in powers of p. By including in the 



previous formulas terms up to the second order in momenta we get [29]: 



nftfr) =n£(o) + flC(p) (2.H4) 

la 
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where 



2„2 



n OT = ^f^^S°^, (2.115) 



,2„2 



and 



s '^ = + Mv^) + v) , (2.117) 



(V-») 2 - F'py-p 
= f + b^V* ^— ^ + (V — > V) . (2.118) 



The coefficients are given by 



^31 _^ l n 2) (2.119) 



108vr V 3 

6 = ( 1 - -ln2 ) (2.120) 

108vr V 3 J y ' 

h = ^ (2-121) 
k 2 = - — (3 + 41n2) . (2.122) 

Since one has, including the one- loop 0(g 2 ) correction in the gluonic la- 
grangian, 

£g = - \ G 2 - - \ G 2 + \ A av ^A hll , (2.123) 

one can compute the one-loop contributions to the Debye and Meissner mass 
from n°°(0) and 11^(0) respectively, according to ( 2.100 ); to the results of 



( |2.115| ) for the diagram of Fig. 2.1a, one must add the contribution ( p. 107] ); 
the two masses are identical for all the eight gluons in the CFL model and are 
given by: 

2 2 

m D = Sri ( 21 " 81n2 ) = 9 2 F 2 , (2.124) 



and 



7T 2 V 36 27 2 



m 2 M = ^t ~^-^ + ^ = (2-125) 



where the first two terms are the result of the diagram of Fig. 2. la and the last 
one is the result of the diagram of Fig. 2.1b (this contribution is called the bare 
Meissner mass in These results have been obtained by the method of the 

effective lagrangian in |2{J and agree with the findings of other authors [45|, 



]46| ]. In particular m 2 M = v 2 mfj, where v, as for the NGB velocity, is equal to 

The Debye and Meissner masses do not exhaust the analysis of the dispersion 
laws for the gluons in the medium. In particular we wish to show the existence of 
a light excitation for small momenta in the CFL phase. To obtain the dispersion 
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laws for the gluons in the small momenta limit we proceed as follows. From 
fl2.112 ) one obtains, for small momenta: 



p°A a-(p-A) . 



(2.126) 



Substituting this result in ( 2 . 1 11| ) taken for v = one gets the dispersion law 
for the time-like gluons A^: 



with 



p° = ±E Ao , E Ao 



a\ = - - 9 - [a — o) 



V3 



\p\ 2 + 



D 



1 + ai 



6A 2 vr v ; 216A 2 tt 2 
The rest mass of Aq is therefore given by 



R 
A 



m D 



y/3(l + ai) 

Since for large values of /i one has g 2 fi 2 3> A 2 it follows that 

m D 



(2.127) 



(2.128) 



(2.129) 



< 



21 — 8 In 2 
6 — A « 1.70 A 



(2.130) 



21 + 16 In 2 

We also define the so-called effective mass (to be denoted m* in the sequel) by 

_ dE p 



dp m* (p) 
in the limit p — > 0. We get, if m* = m*(0), 

y/3 run 



(2.131) 



m A 



3m% w 5.10 A 



(2.132) 



Let us now turn to the dispersion law for the longitudinal and transverse 
gluons; they are obtained by considering (2.111) for u = i and use again (2.126). 
One gets the dispersion laws for the results for the longitudinal and transverse 
cases as follows: 



Ea t 



7T 



+ 



Vi 2 - 1+a2 1 



1 + ai 1 + a-y 



\P 



+ 



D 



1 + «3 

1 + ai ' 3(1 + ai) 



(2.133) 



with 



2„2 



02 



V 9. 



30A 2 tt 2 



(a - 96) 



V 2 9s 



3240A 2 vr 2 



59 



In 2 
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a 3 



A 2 

30A 2 ^ 2 



(a + b) 



3240A 2 vr 2 



41 



112 



In 2 



We get easily the following results 



m D 



73(1+00 



= m R 



(2.134) 



(2.135) 



The results ( 2.130| ) and ( 2.135 ) show that the rest mass is of of the order of 
A and one could therefore wonder if these results are significant, since they 
have been obtained in the limit |p/A| < 1. To estimate the validity of this 
approximation we use the exact result, which can be obtained by equations 
( |2.113| ). Since Il3(p0)0) = 0, the rest mass of the three species Aq, Al, At is 
given by 

m 2 + n 2 (m 2 ,0) = , (2.136) 



with 



n 2 (m 2 ,o) 



+ 



2vr 2 L 
e(x) + e(x 



1 + 



dx 



e(x)[4e(x) 2 - m 2 /A 2 ] 



x + e(x)e(x) 



9e(x)e(x) (e(x) + e(x)) 2 — m 2 /A 2 



with 



(2.137) 



e(x) 
e(x) 

The numerical result of fl2.136| ) is 



m = mR 



1.36 A . 



(2.138) 



(2.139) 



A comparison with ( |2.130| ) shows that the relative error between the two pro- 
cedures is of the order of 18% and this is also the estimated difference for the 
dispersion law at small p. We notice that the gradient expansion approximation 



tends to overestimate the correct result |p3 |. 



2.4.3 Debye and Meissner mass and gluon approximate disper- 
sion laws: 2SC case 

The analysis of the dispersion laws for the glue in the case of the 2SC model is 
complicated by the presence of an unbroken color subgroup. 
The fermionic current appearing in 



A = igAlJH 



can be derived from (2.45). One finds 
A] ( i VhKaclB 







(2.140) 



v A,B=0 ~ V- AaB / 
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where the coefficients Kao.b are given by: 

1 



K AaB = -Tr{\ A X a X B } 



(2.142) 



In this way one can compute the self energy contribution to H^Y . These results 



are contained in [31] and we briefly summarize them. 
We write 

n£(o) + aiC(p). 



with 



and 



K» - 

Km 



2 2 

V 9_ \ ^ v o )( uk 

^2 Z^^ab ' 



TT 2 



(2.143) 
(2.144) 

(2.145) 



It is useful to distinguish gluons of different colours, i.e. the gluons A® (a = 
1, 2, 3) of the unbroken color subgroup SU(2) C and the gluons A® with a = 4, ...7 
and a = 8. 



Colors=l,2,3 

In this case one finds from the self energy diagram of Fig. 2.1 (a, 6 = 1,2, 3): 

JiefVOVV • IV ■ (l + p) + (V <-> V) 



KM 



al) 



+ A 



2vr 3 

Di(£ + p)Di(£) 



D 1 {t + p)D l (£) 



+ 



(2.146) 



As we noted in the previous sections working out the CFL case the extraction 
of the dispersion law in the general case is rather difficult. We consider here 
only the small momentum limit. Working in this limit one gets for the various 
components of the tensor Ii^ u : 



K» = KU0) + 5KUp)= «Op)= 5 ab 



I -»|2 



18vr 2 A 2 



l/ T l 



n£(p) = nS(o) + <m^(p)= 5 ab 5 



rkl 



hi , 



3tt 2 



i + 



6A 2 



and 



KUp) = suZ(p) = s, 



Ok i 



ab 



2 2 

My o fe 
IS^A 2 "^ 



(2.147) 
(2.148) 
(2.149) 



These results agree with the outcomes of |)4|] and [55] . 

We note that there is no contribution to the Debye mass in this case; more- 
over the contribution to the Meissner mass from this term exactly cancels the 
contribution from the tadpole graph. This results reflects the fact that in the 
2SC model the SU(2) C color subgroup generated by the first three generators 
T c (c = 1,2,3) remains unbroken. 
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One can now compute the dispersion laws for the unbroken gluons. From 
the previous formulae one gets: 

£ = -\Fri%, + \v%A%At . (2.150) 



Introducing the fields Ef = F& and Bf = ie ijk F? k , and using (|2~147|) , ( gig ) 



and (2.149) these results can be written as follows (we neglect for the time being 



the 3 gluons and 4 gluons vertices to be discussed below): 

C = \(EtEl - BfBf) + k -ElEt , (2.151) 

with 

* = i&- < 2 ' 152 > 



These results have been first obtained in |35|. As discussed in this paper, these 
results imply that the medium has a very high dielectric constant e = k + 1 and 
a magnetic permeability A = 1. The gluon speed in this medium is now 

oc — (2.153) 



'e\ gn 



and in the high density limit it tends to zero. The one loop lagrangian ( p. 151 ) 
assumes the gauge invariant expression 



C = —FfFf^ (j = 1,2,3) , (2.154) 



provided the following rescaling is used 



A3 _, A3' — 1.3/4 A3 

M - M' = ^4 > 
%o -> x'o = k~ 1/2 x , 

g^g' = k- l / A g . 



Colors=4-7 

For a, b = 4, 7 the polarization tensor is 



KbiP) = iS ab ^^2jd 2 £(v»V»V-£V-(£ + p)+V~v) 



V 

+ n To ,1^ m ) • ( 2 - 155 ) 



0i(* + p)Z> 2 (*) D 2 (£ + p)D 1 {£) 
In the small momentum limit one gets in this case 



KKp) = n» + m» = 5 ab ^L (i + Pb + lg (2 156) 
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DJSOO = SU%(p) = ta b £^P°P l , (2.157) 

and 

n>) = n« ( o ) + ffi >) = { jg(a« + |M + ^_W). 

(2.158) 

We can observe a contribution to both the Debye and Meissner mass in this 
case. 



Color=8 

For the gluon 8, one has 



n 



■ gV 

67T 3 



^2 J d 2 £ (y^vv ■£ V ■ (£ + p) + V <-> f ) 



1 



+ 



Z?i(£ + p) Di(£) D 2 (^ + p)D 2 (^) 
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and, in the small momentum limit, 



(2.159) 



C(p) = n^(0) + OT^(p) 



T oo 



Oil , 



and 



nS(p) = snM 



n£(p) = n&(o) + dna(p) 



vr 2 \ 18A 2 



o i 

54^ 2 A 2 P P ' 

4§ij + 5 ^ + 2 ^ 



(2.160) 

(2.161) 
(2.162) 



Also in this case there are contributions to the Debye and Meissner masses. 
The results are summarized in the following table. 



a 


n°» 


-IL ij (0) 


1 - 3 








4-7 


fm 2 


1 2 


8 


3m 2 g 


1 2 



Table 2.1: Debye and Meissner masses in the 2SC phase; m 2 , = 



These results are in agreement with a calculation performed by [54] with a 
different method. 
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Dispersion laws and light plasmons for gluons 4-8 

To get the dispersion law for the broken color gluons we follow the procedure 
used for the CFL case, Eq. ( |2.108j ) and its sequel. 

Expanding ( 2.112| ) at first order in p we find for b = (3 = 4, 5, 6, 7 



P °A$ = l -p-A? , (2.163) 

while for 6 = 8 

p°A 8 = -p-A 8 . (2.164) 
9 



For colors 4-7 we use eqs. ( 2.156| ) - Q2.158] ) and we get the following dispersion 



laws for the longitudinal and transverse modes: 



-^\p\ 2 + 2A 2 , (2.165) 



Et = ||p| 2 + 2A2 , (2.166) 

The rest mass for these gluons, in the gradient expansion approximation, is 
given by 

™A L = m A T = v^A. (2.167) 

To estimate the validity of this approximation we can use the exact result, 
obtained analogously to the method described for the CFL case. The rest mass 
of the three species Aq, Al, At is given by 

m 2 + n 2 (m 2 ,0) = , (2.168) 

with 



(x + v x 2 + l) 2 — m 2 v x 2 + 1 



(2.169) 

The numerical result of ( |2. 168 ) is 



m = mR = 0.894 A . (2.170) 

A comparison with ( |2.167| ) shows that the difference is of the order of 40 — 50% 
and this is also the estimated difference for the dispersion law at small \p\. 
We notice that also in this case, as for three flavors, the gradient expansion 
approximation tends to overestimate the correct result. 

In the case of the gluon of color 8 we work as for the colors 4-7. We obtain 
the following dispersion laws for the longitudinal and transverse modes: 



El = \/ 270^bT + ^, (2.171) 



1 m 2 . 



Et = \/-3o^l#| 2 + ™ 2 M) (2-172) 



43 



where the Meissner mass is in table 2.1. From these equations we see that 



m M ■ (2.173) 



The large rest mass of the 8th gluon is an exact result, independent of the small 
momentum approximation and is due to the fact that, in this case 

2 2 

n 2 (m 2 , 0) = - = -m 2 M . (2.174) 

2.5 One-loop corrections to three- and four-point 
gluon vertices in the 2SC model 

To complete the effective lagrangian for the eight gluons in the two-flavor model 
one can compute the one loop corrections to the gluon vertices (3 gluons) 
and 1^4 (4 gluons). We shall not consider in the sequel a possible light glueball 
which is discussed in [56]. Therefore we write the full gluon lagrangian as 

where and C 1 ^ are the one loop lagrangian terms for the three and four 
point gluon vertices respectively. For the three point gluon vertex only the 
diagram in Fig. 2.2 survives at leading order in 1///; for the four point gluon 
vertex only the diagram in Fig. 2.3 survives in the same limit. 



Fig. 2.2 Three gluon vertex. Dotted lines represent gluon fields; full lines are 
fermion propagators. 

To start with, let us consider the three-point function. At the tree level (£ = C ) 
the contribution to the lagrangian can be written in the form: 

4) = -gfabcA^d^ . (2.176) 

At one loop (C = C 1 ) we have to distinguish between the contribution of the 
diagram involving the gluons in the unbroken SU (2) gauge group, which we call 
£L and those involving gluons corresponding to broken generators, which we 
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Fig. 2.3 Four gluon vertex. Dotted lines represent gluon fields; full lines are 
fermion propagators. 



call £^ 2 - So the one- loop correction at the three gluon vertex may be written 
as follows: 

£(3) = ^(3),1 + £(3),2* (2.177) 

For the SU(2) contribution our result is as follows: 

£(3),i = -gkfabcA^Ald^A^d^i + SfuSrf] , (2.178) 
g 2 up 

with a,b,c 6 {1,2,3} and A; = 0-1-9 • This term is consistent with the fact 

that color superconductivity leaves invariant the subgroup SU(2) C generated 
by the SU(3) C generators T a , a=l,2,3. It can be also obtained in a simpler 
manner, by requiring gauge invariance for the SU (2) gluons. We do not write 
down explicitly £^ 2 whose expression can be found in [j3l]| . 

The same results are found for the four gluon vertex; again the structure of 
the vertex coincides with the one dictated by gauge invariance for the gluons 
1,2 and 3, while this invariance is broken for the other generators. 



2.6 Gap equation in HDET 

Let us conclude this section by deriving the gap equation in the High Density 
Effective Theory. Let us begin with the CFL case; to get the gap equation one 
can use the effective action formalism of Cornwall- Jackiw-Tomboulis (5?]]. In 
this context it corresponds to the calculation of the grand potential 



n 14 ^ 



2 {2ir) 4 



J d 2 l { Tr In [s^ ee {e)S(£)] - TrS{t)Z{£)} . (2.179) 



Here 5 1 = and S, ee = S A g j ree are respectively the inverse fermion 
propagator and the free inverse fermion propagator and X = S" 1 — Sj ree is 
the proper self-energy. By extremization of the density energy with respect 
to the gap one gets the Schwinger-Dyson equations; the fermion equation is 
the desired gap equation. In our case the truncated Schwinger-Dyson equation 
assumes the form 

S^p = 0) - S^jrJp = 0) = 
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A-rr^g 2 



1 U ab p. 



aAC 



Here 



J 



aAB 



9 



i V^flAaB 




iS C D(£)J: 



o 



bDB 



with liAaB given in ( 2.103 ). We assume a fictitious gluon propagator 



A 2 



(2.180) 
(2.181) 

(2.182) 



which corresponds to a local four fermion interactions similar to the Nambu- 
Jona Lasinio model H. One obtains two equations: 



2„2 



A 

A 9 



Zip 9 

4A 2 7T 3 

2ifj 2 g 2 A 
3A% 3 



^I(A 9 ) - p ( A) 



where 



1(A) 



I 



d 2 i 



V -IV ■ 



A 2 + it 



i 2-7T arcsinh ( - — - 
A 



(2.183) 



(2.184) 



We fix A/g by using the same equation and the same model for the gluon 
propagator at fx = 0, where the effect of the order parameter, in this case the 
chiral condensate, is to produce a constituent mass M for the light quark. This 
assumption gives the equation 



1 



p 



3A 2 vr 2 Jo + M 2 ■ 



(2.185) 



For M = 400 MeV and the cutoff K = 800 MeV (700 MeV) we get A/g = 181 
MeV (154 MeV). The corresponding values for the gap parameters are given in 
tables 2.2 and 2.3. 



K = 5 + n 


A 


A 9 


A 9 /A 


800 MeV 
700 MeV 


66 MeV 
77 MeV 


-150 MeV 
-189 MeV 


-2.27 
-2.45 



Table 2.2 A, A 9 and A 9 /A in the CFL model; fi = 500 MeV. 



K = 6 + n 


A 


A 9 


A 9 /A 


800 MeV 
700 MeV 


25 MeV 
51 MeV 


-57 MeV 
-114 MeV 


-2.28 
-2.24 



Table 2.3 Same as in table 2.2, with \x = 400 MeV. 

In order to obtain these results we have used the same cutoff K for the two 
gap equations at fj, = and /* / 0; this fixes the value of 5, the cutoff on the 
residual momentum along the Fermi velocity at fi ^ (|/||| < 5): 

S = K-fi. (2.186) 
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A similar approach has been used by other authors, see for example the review 
1 14] and references therein. Their calculation differs from the present one for 
two reasons: they include the 0(1/ fx) corrections to the gap equations and use 
a smooth cutoff 

with A = 800 MeV and various values of v. They obtain at /x = 400 MeV 
A = 80 MeV and A 9 = -176 MeV while for p = 500 MeV A = 109 MeV and 
Ag = —249 MeV. Note that in any case the ratio Ag/A ~ —2 which means 
that the CFL ansatz is indeed justified by this analysis. 



The analogous of eq. (2.183) for the 2SC case is 



that can be solved explicitly: 

IAI = - 



3A^ arCSmh I IAI 



sinh 



3AV 



(2.188) 



(2.189) 



With the same values of A/g (A/g = 181 MeV and 154 MeV) we find the values 
for the gap parameter of the 2SC model reported in tables 2.4 and 2.5. 



K = S + fi 


A 


800 MeV 
700 MeV 


88 MeV 
105 MeV 



Table 2.4 Value of the parameter A in the 2SC model; fi = 500 MeV. 



K 


A 


800 MeV 
700 MeV 


39 MeV 
68 MeV 



Table 2.5 Same as in Table 2.4 with \x = 400 MeV. 

2.7 Appendix. A few useful integrals 

We list a few 2-D integrals used to obtain the results of this Section. Let us 
define 



I 



d N e i (-i) n+1 nf T(n + 1 - f ) 

(y ■ £ V ■ £ - A 2 + ie) n+l ~ n\ &2n+2-N ■ 

therefore, for N = 2 — e and with 7 the Euler-Mascheroni constant, 



(2.190) 



2Z7T o 

1q = M7rln7rA + m^f , 

iir in iir 

1 = + A2 ' 2 = ~2A4 ' 3 = + 3A6 ' 



(2.191) 
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Moreover we define 
/, 

We get 

h,i = 

h,i = 

h,i = 

h,2 = 



d 2 



{V • I V ■ t - A 2 + ie) n (V ■ t V ■ I - A' 2 + ie) 
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The following integral is used in the context of the 2SC model: 

(V ■ if d 2 i 



2.192) 



h 



(V -IV -if 



(2.193) 



The divergence is cured going to imaginary frequencies and finite temperature 
(£o = iu> n , uj ti = irT(2n + 1); T — ► 0, fj, — ► oo): 



Ii = 2nTi 



Other divergent integrals, such as 

i - 



2nTi ( — — tanh 

T 2T 



d 2 £ 



-2m . (2.194) 



(V ■ IV ■ £)(V ■ IV ■ I - A 2 ) 



(2.195) 



are treated in a similar way. 

Finally, useful angular integrations are: 



dv , , V k 
— V J V K = — 

4vr 3 



An 15 



(2.196) 
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Chapter 3 

Theoretical developments 



3.1 Advances in the CFL model 

In the last three years there has been a tremendous activity in the study of 
QCD at high density. To review these results is beyond the scope of this paper 
and I refer to some excellent review papers that have been already published 



1 14]. I will survey here only some of the most important developments that 
in a way or in another are related to our subject, i.e. the effective lagrangian 
approach. To begin with, in this paragraph I report on some of these advances 
in the Color-Flavor-Locking model. 

3.1.1 Would-be NGBs masses 



In par. L5 the effective theory for the NGBs in the CFL was presented^. In the 



CFL model the Nambu Goldstone Bosons are the only elementary excitations 
that are massless in the limit of zero Dirac quark masses. To have a more 
realistic description one should however include the effect of Dirac quark masses, 
which implies that also the octet of would-be NGBs acquire mass (the superfluid 
mode remains massless). The Dirac masses induce a symmetry breaking term 
in the lagrangian that becomes: 

F 2 / \ 
C eff = — Tr [dtEdttf - v 2 |VS| 2 J - c [detM Tr (M _1 S) + h.c] (3.1) 



where M = diag(m u , m^, m s ). The value of the parameter c is J45|, [47|, |4^, [49 



5C] : 

3A 2 , , 

c~. (3.2) 

The masses of the different bosons in the octet depend linearly on the Dirac 
quark masses and have a pattern different from the analogous octet of pseudo 
NGBs of the chiral symmetry; for example one finds: 

m 2 ± = -r^m s (m u + m d ) . m 2 K± = ^m d (m u + m s ) . (3.3) 



In the 2SC model there is one NGB and three gluons remain massless. The effective 
theory for the 2SC model has been developed in |H| . Its coupling to electromagnetic field and 
weak currents is in |B8|. 
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Note that 

4± * , (3.4) 

i.e. the kaon is lighter than the pion. 

Numerically these masses are much smaller than the corresponding masses 
at zero densities, due to the suppression factor A///; for example the kaon mass, 
for fi ~ 500 MeV is in the range of ~ 10 MeV. 

3.1.2 Electromagnetic interactions 

The diquark condensate breaks U(l) em , which is clear as it couples a quark pair 
with a non zero total charge. In presence of spontaneous breaking of electro- 
magnetism photons acquire a Meissner mass and ordinary conductors become 
superconductors. In the case of color superconductors, however, the situation 
is different. Here a new U(l) subgroup appears, arising from a combination of 
color and flavor; it remains unbroken and takes the role of electromagnetism in 
the high density phase, while the associated vector boson remains massless j| . 
The situation is similar to the standard model of the electroweak interactions, 
where the photon arises as a combination of the two neutral bosons related to 
two diagonal generators of the SU(2) x U(l) electroweak group. 

To determine the unbroken U(l) group, let us show that the condensate 

rr[V> T e/Vei] , (3.5) 
see Eq. ( [2. 33 ), is invariant under the U(l) transformation 



with 



ipai -> U^lppj (3.6) 



Uff = (exp { ie [1 ® Q + T ® , (3.7) 



where Q = diag(2/3, —1/3, —1/3) is the quark charge matrix in the flavor space 
and T a p = l/3(— 2, 1, 1) acts on color indices. In infinitesimal terms: 

^ _, ip + , (3.8) 

where 

5ip ai = i e [ip aj Qij + T al3 ippj } , (3.9) 
To prove the invariance under the U{\) group generated by 

Q = 1®Q + T®1 (3.10) 

it is sufficient to note that 

5Tr[^j T e I i/je I ] = 2ie{[Tr{ipTip + tpi/;Q) - Tr(i/))Tr(Tip + ?f)Q)]} = 0, (3.11) 

where we have used (|2.35| ). Therefore the unbroken electromagnetism arises 
from the combination of the original U(l) em with the U(l) subgroup generated 
by a color hypercharge. 
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Let us now consider the vector boson taking the place of the photon. It 
would appear in the lagrangian through the covariant derivative 



Dp = a M - i [e Q + g G*VST) (3.12) 
(T 8 = V3T). Writing 

An = cos 9A fl - sin 9G\\ , (3.13) 

G 8 = + sin 9A fl + cos 9G 8 ^ , (3.14) 

one can see that the field coupled with a strength proportional to Q is A^, 
obtained from these equations with 

cos9 = ^J^=. (3.15) 

Therefore the field 

A^ = cos 9A^ + sin 9G* (3.16) 
takes the role of the massless photon field, while 

G 8 = — sin 9 A^ + cos #G 8 (3.1 7) 

is a massive gluon. We note that, since </>e, the angle 9 is small, 9 ~ and the 
effective photon is mostly made by the ordinary photon. We also note that the 
electromagnetic charges Q of the quarks are produced by combining two sets of 
charges, the charges <7j = (2/3, —1/3, —1/3) and the charges t a = 1/3 (—2, 1, 1); 
the quarks have therefore always integer electromagnetic charges, as in the 
old Han-Nambu model. The unit of the electric charge, corresponding to the 
absolute value of the electron charge, isf] 

69 (3.18) 



2 



9 2 + % 

The quarks up with colors r, g, b have respectively charges 0, +e, +e; the quarks 
down (or strange) with colors 1,2,3 have respectively charges — e, 0, 0. It is 
worth noticing that also the gluons have integer electric charge Q, exactly as 
the octet of vector bosons of low density QCD. The discussion of the electroweak 



couplings in the context of the effective theories can be found in [58|. 



3.1.3 Quark-hadron continuity 

Let us summarize once again the properties of the low energy spectrum in the 
CFL model. Besides the octet+singlet quasi-particles of spin ^ with integer 
charge there is an octet of massive vector states, again with integer charges and 
an octet of pseudoscalar NGBs: they are massless in the limit of zero quark 

2 The electron charge is less than e because the electron couples to the photon with a 
strength reduced by cosd. 
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masses, but acquire masses if the chiral symmetry is explicitly broken; finally 
there is a massless NGB associated with the breaking of U(1)b that remains 
massless even in presence of quark masses. There is also a light pseudoscalar 
bosons that is never massless: it is associated to the symmetry U(1)a which 
is broken by instantons. If we now compare these results with the symmetries 



of the 3-flavor hypernuclear matter [59] we see a considerable similarity. On 



one side we have the integer-charge low lying baryonic octet N, S, A whose 
members can pair to produce SU(3) condensates, e.g. AA, EE etc. breaking 
the baryonic number; on the other side there are the octet+singlet quasiparticles 
of spin 1/2 with integer charges as well. On the nuclear side we have the low 
mass pseudoscalar octet of pions, kaons etc. and on the other side the octet 
of NGB of the broken symmetry of the CFL phase. Finally to the low mass 
J p = 1~ octet of vector resonances (p, K* , etc.) would correspond the massive 
gluons. All that means that probably the two phases are contiguous and one 
passes from one to the other without phase transition. A further example will 



be discussed in paragraph 3.2.2 



3.1.4 Asymptotic analysis 



We have derived in the paragraph 2.6 the gap equation in the Nambu-Jona 



Lasinio model where the gluon propagator is substituted by a local four-fermion 
coupling. Eq. ( |2.183; ) has the approximate solution in the CFL model 

f 3A 2 7r 2 1 

A«25exp|-^p-J . (3.19) 
A more refined treatment would involve the use of the gluon propagator and 



has been performed in |p0|, |pl]1 (see also [S2| and Barrois and Frautschi in |TJ). 



The result of this analysis is given by the formula: 



H ___ f 37T 2 \ 

V2g) ' 



A «r_exp . 



One should note the different behavior 

A ~ exp(— const / g) (3.21) 

instead of A ~ exp(— const /g 2 ) between (|3.20| ) and (|3.19| ), which is due to the 
fact that in the complete treatment the gluon propagator at high density has a 
collinear divergence at zero scattering angle. We can qualitatively understand 
the behavior ( |3.21| ) as follows [ 14 ] . The collinear divergence should be regulated 



by the gluon rest mass that, as we know, is of the order A, see par. |2.4,2 . 
Therefore the gap equation has the form 

A oc 5 V / dlodO ao * A2 . (3.22) 



+ a 2 v?e + A 2 



In this way one gets the relation 



lc<3 2 (ln^ (3.23) 



and therefore eq. ( p. 21 ). 
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3.2 Other models 



3.2.1 One flavor 

In the real world m u ~ ~ might be a good approximation; therefore, if 

A > m 2 J(2/i) , (3.24) 

the color flavor locking could be a good approximation as well. The reason for 
the condition ( |3.24| ) can be understood as follows (see |3^] and the discussion 
on the strange quark mass in j| and (7j). The main effect of the strange quark 
mass is kinematic and produces a shifting between the Fermi momenta of the 
u, d quarks, pf — M, and that of the strange quark, pf = y^ 2 — m 2 ~ /i — ^-j. 
The difference between the two momenta must be small for the pairing to occur, 
in particular it must be smaller than the smearing of the Fermi surface, of order 
A, which explains ( |3.24| ). 

On the other hand, if m s > ^2[iA, one has to consider the role played 
by the strange quark explicitly. Quantum Cromo Dynamics at large chemical 
potential and small T with one flavor (strangeness) has been considered in 
]43[ ]. Also in this case the theory exhibits color superconductivity, however 
the symmetry breaking condensate, being antisymmetric in color, must be in 
an angular momentum 1 state, because there is no other way to obtain an 
antisymmetric wavefunction for the strange quark pair. In the non relativistic 
limit the computations are relatively simple fj|; there are two possibilities for the 
order parameter A?, i.e. 

<^ T CVA^> , ^ T Cf\ a A ^); (3.25) 

here (a = 1, 2, 3) are the three Gell-Mann matrices antisymmetric in color 
(i.e. A2, A5, A7); the two terms in ( 3.25| ) are distinct because in the non rel- 



ativistic limit spin and orbital angular momentum are separately conserved. 
Four possible cases were considered 



Af 


= A5f, 




A- 


= A5 a3 5 i3 




A? 


= A5 a3 (6 n 


+ ik2) , 


A? 


= A(5 al 5 ll 


+ 5 a2 8 i2 ) 



(3.26) 

called respectively Color Spin Locking (CSL), 2SC, A and planar phase. For 
the former condensate in ( 3.25| )^|the results are as follows. First one determines 



the value of the parameter A solving the Schwinger-Dyson equation for the four 
cases in ( |3.26| ); one obtains that the largest gap corresponds to the CSL case. 
Next one looks at the grand potential and observes that the solution minimizing 
it corresponds again to the Color Spin Locking phase. One can conclude that 
the CSL phase is indeed favored for the case of one- flavor QCD at high density. 

3 The analysis of the ultrarelativistic limit does not change the conclusions obtained by the 
non relativistic approximation. 

4 The latter condensate produces either smaller or anisotropic gaps and the analysis is more 
involved, see 1431 for further discussion. 
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The global symmetries of the CSL phase agree with the expectations based on 
Nf=l QCD at high density, and therefore one can argue that also in this case, 
as for CFL model with Nt = 3, low and high density phases are continuously 
connected. 



3.2.2 QCD at large isospin density 

It would be desirable to confront the analytical results obtained by different 
approximations for QCD at high density with the lattice QCD calculations. 
However this comparison is at present impossible because current methods to 
evaluate the QCD partition function use path integral with a complex fermion 
determinant. In the quenched QCD approximation one ignores the fermionic 
determinant, but this procedure leads to inconsistent results in the case of 
finite baryonic chemical potential fig |p3|, [S4|. In absence of a consistent 



approximation scheme, one may hope to get a clue to the behaviour of quark 
matter at finite \xb by considering QCD at finite fii, where fij is the chemical 



potential associated to the third component of isospin, 1%. In |65] the case of 
finite fii and [is = was considered, which, though unrealistic, can be of some 
utility because it should be accessible to lattice calculations [pq]. The results 



found in 6E] are as follows. When < the vacuum does not differ from 
the [ii = case, because to excite a ir costs — If is small, but 

larger than m n , the excitation of the pions is favored and the ground state is 
a pion superfluid (with one massless pion and the other two massive). This 
state is characterized by the order parameter < 7r~ >^ 0. For large one 
goes to a region characterized by the order parameter < u'j^d >^ 0. Since the 
two order parameters < 7r~ and < u^d have the same quantum 
numbers and break the same symmetries it can be conjectured that there is no 
phase transition along the fij axis from the small to the high chemical potential 
[|65f and one can move continuously from the hadron to the quark phase. This 
hypothesis is similar to the conjecture of quark hadron continuity discussed in 



the CFL model, see paragraph 3.1.3 



3.3 Transport theory for quarks in dense matter 

The transport theory for colored particles in a Yang Mills external field can be 
used to derive quasi-particle properties for QCD at high density. This develop- 
ment not only allows a different derivation of known results, but also sheds a 
new light on the dynamics of the quasi-particles. 

The classical equations describing the motion of a particle in colored external 
field have been derived by WongP] f67|| . Here we derive the Wong equations by 
the Fok Schwinger proper time method that presents the advantage of being 
covariant. 

To begin with, we write the Dirac equation for the Feynman propagator 

(ilftx — iti)Sf(x, x') = 5 4 (x — x') 

5 The original Wong derivation used Dirac equation. The quadratic form of Dirac equation 
was used by Heinz [|68f to disentangle quark spin. 
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in an operatorial way 

(ity — m)Sp = 1 . 

We have 

g 1 itf + m —2m il} + m 

F iljt — m —2m (il?) 2 — m 2 —2m ' 

with 

_ -2m l_ 

(j^)2 _ m 2 = ■ 

Let us interpret H, given by 

_ (iff) 2 -m 2 
* " -2m ' (3 ' 27) 

as a hamiltonian; with r proper time of the quark we have: 

HV = i — , (3.28) 

i.e. 

dU 

i-g F = EV, 

where 

U = U(t) = e~ iHr 
is the propagator. In the coordinate representation one has 

U(x,x',t) = < x\U(t)\x' > . 

Since one has 

1 f o 
G = — = - i / dTe~ lHr , (3.29) 

we also have 

id T <x\U(t)\x' > = < x\HU\x' > =< x\UU ] HU\x' >= 
= <x(t)\H(x(t),tt(t),Q(t))\x'(0) > , 

where 

^ = iD» = + gA^ (3.30) 
is the kinetic momentum, and 

\x( T ) > = U\t)\x(0) > , 

\x> = \x(0) > . 

Moreover 

U\t)HU{t) = H(x(t),7t(t),Q(t)) 
is the hamiltonian and 



7t(t) 
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x(t) = U\t)xU{t) , 

Q a (r) = U\r)Q a U{r) 

are operators in the Heisenberg representation. 

The Heisenberg equations of motion are as follows: 

If = ilH - *"> ■ 



dlt^ 

~d~7 

dQg 

dr 



i[H, ^] , 
i[H, Q a ] ■ 



(3.31) 



Since the following expression holds for H: 



H 



-1 



— (f-m){if + m) = 

[p^ + gA^Y + g -a-F 



-1 
2^i 



m 



(3.32) 



the Heisenberg equations are transformed as follows. The first equation in 
( 3.32| ) becomes 

dx* ir^ 



dr m 

as expected. The analogous of the Lorentz force is as follows: 



(3.33) 



dir» 

~dV 



— ( tt x F^ + F A/ %a - -D^tr ■ F) 



2m 



1 



while the dynamics of the color charge is governed by the equation 



dQ a 
dr 



2m 



1 



Since in the classical limit all the operators commute, one has 

[n,F] = 0, 



[Q a ,Q h 



ifrfabcQ C -»■ 0. 



(3.34) 



The spin terms can also be neglected and one obtains the Wong equations 
(FV = -QapaXy, and Q a analogous of -A/2 in the classical limit p§): 



m 



m ■ 



dx^ 
~dV 
dir" 
~dT 

dQg 

dr 



7T' 



g Q a F a » x 



7TA , 



-gfabcTTX A bX Q C . 



(3.35) 
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Let us now consider the classical transport equation for the one particle 
distribution function f(x,p, Q) which is the probability for finding the particle 
in the state (x,p,Q): 



df(x,p,Q) r , 
m — = C[f](x,p,Q) , 

where C[f] denotes the collision integral. Putting 

Dp = dfj, — gf abc Q c Affi , 



(3.36) 



(3.37) 



from equations ( |3.35| ) in the approximation C[f] = we get therefore the 
Boltzmann equation in the collisionless case: 



) 



Notice that in (3.38|) the velocity 



de(p) 



(3.38) 



(3.39) 



is determined by the dispersion law of the quasiparticle. 

Eq. ( [3.38 ), together with the Yang Mills equations, that we do not write 
here, represent the complete set of Boltzmann- Vlasov equations for the non 
abelian case and have been considered in [69] to derive the hard thermal loops of 
QCD. In [70] Litim and Manuel apply it to the two-flavor color superconductor. 
The method can be summarized as follows. First of all one writes the colored 
current associated to the motion of the quasi-particle: 

W = sE J j0tf d Qv»Q a f{x,P,Q) , (3.40) 



where the sum is over the helicities and 

1 



[l,v). Let us call 



/ 



eq 



the Fermi equilibrium distribution function and let us assume a small, pertur- 
bative deviation from equilibrium: 

f = f e9 + gSf . 

One also defines the color density 



J a (x,p)=g J dQQ a f(x,p,Q) , 

so that ( EOH) 

can be written as follows: 

d^p 



(3.41) 



{2irH) 



rV M J a (x,p) 
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Expanding (|3.38|) in g one obtains the transport equation for the color density: 



D t +v-D 



J(x,p) = g 2 N f v-E 



df 



eq 



de 



The formal solution of ( 3.42 ) is 



K = 9 2 N f 



d 3 p dry 



4„, . Jfeq 

v,G ab v.E b (y) 



3 



de 



(3.42) 



(3.43) 



where G ab =< y\(v ■ D)~ l \x > a b is the Green function. We observe that, since 
D depends on the gluon fields A, is a functional of A. This treatment is 
analogous to the derivation of the Hard Thermal Loop and the Hard Dense 



Loop approximations from the transport theory [69], [71], [72]. For the non 
relativistic and Abelian case it reduces to the transport theory for the BCS 
superconductor [f73|] . 

All the relevant information is contained in the functional J [A] because the 

effective action T can be obtained from J\A] = — -— r and all the diagrams can 

L J SA 

be obtained by deriving T. As a matter of fact, expanding J [A] in powers of 
An one gets 



T-rao a v | -paoc a v a A 



(3.44) 



Therefore from ( |3.43 ) one can obtain the various tensors, in particular the 



polarization tensor n^, that in |7C] is computed for the 2SC model. 



There are several advantages in computing the polarization tensor by this 
method. First of all one obtains at once not only the results at T = 0, but also 
the thermal dependence. For example the T— dependence of the Debye mass is 
given by 

,df eq g 2 n 2 N ft /A T N 



m 2 D 



where 



9 2 N f 
' 2tt 2 



I" I rp 1 



dpp 



de 



1 



2^ 2 



In I - 



(3.45) 



1 , 2 n df 
dpp v 



eq 



de 



M Jo 

One can check that I n (oo, 0) = 0, which implies that there is no Debye screening 
in the 2SC model at T = 0, a result we already knew. For small temperatures 
one gets 

„2,, 2; 



m 2 D 



2vr 2 



gViV/ / 2ttA _ a/t 



(3.46) 



From the real part of H^f u one can get the dispersion laws; the resulting expres- 
sions are involved, but can be expanded for small momenta, similarly to the 
approximations presented in the par. |2.4.2 . On the other hand the imaginary 
part of n^, gives the Landau damping, that, however, is absent for small ener- 
gies, i.e. plasmon and transverse excitations are stable provided EL^ik) > k. 
In |7(]] a comparison with the results obtained by quantum field theory at T ^ 
[54] is given and an agreement is found. The method based on the transport 
equation is however more powerful because it allows the calculation of other 
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significant quantities beyond the leading order, such as, for example, the damp- 
ing rate, or transport properties such as thermal and electrical conductivities or 
shear viscosity that are dominated by the excitations of the NGB or the gapless 
quarks. 



59 



Chapter 4 

The LOFF phase 

4.1 Overview of the crystalline superconductive (LOFF) 
state 

In nature flavor symmetry is broken not only explicitly by quark mass terms, 
but also by weak interactions. Therefore in the applications of the color su- 
perconductivity one has to take into account this symmetry breaking; for ex- 
ample in compact stars, considering only two flavors, isospin is broken by 
Sfj, = \i u — [id 7^ 0, due to the process: 

d — > uev . (4.1) 

When the two quarks in the Cooper pair have different chemical potentials, the 
vacuum is characterized, for certain values of 5/i, by a non vanishing expectation 
value of a quark bilinear breaking translational and rotational invariance. The 
appearance of this condensate is a consequence of the fact that in a given range 
of 5y Iffi l, the formation of a Cooper pair with a total momentum 

p 1+ p 2 = 2q^0 (4.2) 



is energetically favored in comparison to the BCS state characterized by q 
k + (—k) = . For the two flavor case one finds [35| 



Sfi e [<fyii,<5/i 2 ] , 

fyti = 0.71A , Sfi 2 = 0.74A, Sfij < y. (4.3) 

where A is the BCS gap. A similar phenomenon was already observed many 
years ago in the context of the BCS theory for superconducting materials in 
presence of magnetic impurities by Larkin, Ovchinnikov, Fulde and Ferrel and 
the corresponding phase is named LOFF state p0| . 

The exact form of the order parameter (diquark condensate) breaking space- 
time symmetries in the crystalline phase is not yet known. In |35j the following 
ansatz is made: 

A(x) = Ae i2 ^ = Aj 2qH - s . (4.4) 
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Fig. 4.1 LOFF kinematics. The Cooper pair has a total momentum 2q ^ 0. 



The value of |<j| is fixed by the dynamics, while its direction n is spontaneously 



chosen. In the sequel I will assume the ansatz (4.4) as well. The order parameter 



( |P| ) induces a lattice structure given by parallel planes perpendicular to n: 

n-x = — (k = 0, ±1, ±2, ...) . (4.5) 

q 

We can give the following physical picture of the lattice structure of the LOFF 
phase: Due to the interaction with the medium, the Majorana masses of the red 
and green up and down quarks have a periodic modulation in space, reaching 
on subsequent planes maxima and minima. 

The Cooper pair can be formed only if quarks are in the antisymmetric 
channel, where there is a color attractive interaction; therefore it must be in 
antisymmetric flavor state if it is in antisymmetric spin state (S = 0). The 
condensate has therefore the form 

- < Ole^e^^^C^^lO >= 2T^e 2 ^ . (4.6) 

The analysis of [35| shows that, besides the condensate (WM) (scalar condensate), 



another different condensate is possible, i.e. one characterized by total spin 1 
(vector condensate) and by a symmetric flavor state: 

i < Ol^e^sVn^C^V^lO >= 2T L B e 2i ^ . (4.7) 

Note that in the BCS state the quarks forming the Cooper pair have necessarily 
5 = 0. 

It goes without saying that the hypotheses of |35| are rather restrictive, as 
these authors assume only two flavors and make the ansatz of a plane wave 



behavior, Eq. (4.4). In any event their results are as follows. Assuming a 
point-like interaction as the origin of the Fermi surface instability, the LOFF 
state is energetically favored in the above-mentioned small range of values of 5fi 
around 5fi ~ 0.7A. The actual value of the window range compatible with the 
presence of the LOFF state depends on the calculation by which the crystalline 



color state is computed. While the determination of the interval (|4.3|) is based 
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on a local interaction, assuming gluon exchange, as in |36|], the window opens 
up considerably. The value of the gap parameters will be discussed below, at 
the end of par. 4.3. 



The order parameters ( f4,6|) and (4.7) spontaneously break rotational and 
translational symmetries. Associated with this breaking there will be Nambu 
Goldstone Bosons as in a crystal; these quasi-particles are known as phonons. 
Let us discuss them in some detail. 



4.2 Phonons in the LOFF phase 

For a generic lattice structure it is known that there are three phonons associ- 
ated to the breaking of space symmetries. However one can show ]32|] that in 
order to describe the spontaneous breaking of space symmetries induced by the 
condensates fl4.6|) and ( |4.7| ) one NGB is sufficient. The argument (see Fig. 4.2) 
is as follows^. Rotations and translations are not independent transformations, 
because the result of a translation plus a rotation, at least locally, can be made 
equivalent to a pure translation. Let us discuss the consequence of this fact 
and consider the goldstone field associated to the spontaneous breaking of the 
rotational and translational invariance. In the present case the NGB is a long 



Fig. 4.2 In the point P the effect of the rotation n — > ft' and the effect of the 
translation x — > x + a tend to compensate each other. 

wavelength small amplitude variation of the condensate A(x) — > A(x)e' 1 ^^ , 
with 

cf)/f = 2q{n + 8n) ■ (x + 5x) - 2qn ■ x = 2qR ■ x + T - 2qn ■ x , (4.8) 

where we have introduced the auxiliary function T, given by T = 2qR- Sx. 

Now the lattice fluctuations 4>/ f must be small; from this it follows that T 
and R are not independent fields and therefore T must depend functionally on 

1 A similar argument is used in [|74| in a different context. 
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R, i.e. T = F[R], which, using again ( |4.8| ), means that 

$ = 2qn ■ x + <f>/f = 2qR ■ x + F[R] = G[R, x] . (4.9) 

The solution of this functional relation has the form 

R = h[$] (4.10) 

where h is a vector built out of the scalar function $. By this function one can 
onlyQ form the vector ; therefore we get 



R 



IV$| 



(4.11) 



which satisfies \R\ = 1 and < R >o= n. In terms of the phonon field eft the field 
R is given at the first order in <fi by the expression 



R = n + 



2fq 



V</> — n(n ■ V0) 



(4.12) 



In order to derive the form of the interaction between the NGB and the quarks 
let us now write the effective lagrangian in the LOFF phase. 



4.3 HDET approach to the LOFF state 

To take into account the two non vanishing vacuum expectation values we add 
the term: 

1 X e a ^f a (x)C (A(% + a ■ nA^4) ^(x) 



2 - riut\—/- ^ — ~*j ■ — — 13 j 
+ h.c. (4.13) 

to the lagrangian. 

Let us first consider the scalar condensate: 

£a = -^e^e^e^Kx^fsix) -(L^R) + h.c. . (4.14) 

Here ip{x) are positive energy left-handed fermion fields discussed in the pre- 
vious section. We neglect the negative energy states, consistently with the 
assumptions of HDET. 

In order to introduce velocity dependent positive energy fields 



2 In principle there is a second vector, x, on which R could depend linearly, but this pos- 
sibility is excluded because R is a vector field transforming under translations as R(x) — > 
R'(x') = R(x). 
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with flavor i, we have to decompose the fermion momenta: 

Pj=HVj+?j (i = l,2). (4.15) 

Therefore we have : 



£a = — ^exp{ix ■ a(vi, v 2 ,q)}eij 

x e a ^ +t ^ i;ia (x)C^ +t ^. ;jp (x)-(L^R)+h.c. (4.16) 



where 



d(vi, v 2 , q) = 2<f - MiUi - • (4.17) 

We choose the z— axis along q and the vectors vi, v 2 in the x — z plane; if a±, a 2 
are the angles formed by the vectors v±, v 2 with the z— axis we have 

®x = — Ml sin a\ cos 4>\ — \i 2 sin a 2 cos 4> 2 , 
a y = — fjii sin ai sin </>i — fi 2 sin «2 sin cp 2 , 
a z = 2g — /^icosai — ^2Cos«2 • (4-18) 

In the fix, H2 -^00 limit the only non vanishing terms in the sum correspond 
to the condition 

a x = a y = 0. (4.19) 



Let us introduce 



_ Ml + M2 

f/i = (4-20) 



ai = Q2 + vr + ( - ) , (4.21) 



Eq. (p^| ) impl ICS 

«i = a 2 + tt + O ( 

i.e. the two velocities are opposite in this limit. 

The limit /1 — > 00 does not imply a z = 0; however if we might take the 
limit q ~ <5/i — > 00 as well, we would also have a 2 ~ 0; we will make this 
approximation [34], which is justified by the results of [35]: q ~ 1.25 fi and 
e)/i ~ 0.7Abcs 3> A( s ' u ). This approximation makes the loop calculations 
much simpler than the analogous in [35] and therefore it may be particularly 
useful in complex calculations. In any case corrections to this approximation 
can be implemented. We have up to terms O 

<5/i q 2 - <5/i 2 
cos a± = h 



qfj, 

cosa 2 = + ^ + ^!, (4.22) 
q qn 



which, together with ( 4,19| ), has the solution 

5fi e 

7 " 2 

ai = a 2 + tt - e (4.24) 



ot 2 = q = arccos — — - , (4.23) 
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with 



1 1 



(4.25) 



Therefore, as anticipated, e = O ( — 1 and in the limit u — > oo the two velocities 

V/V 

are almost antiparallel: 
Putting 



V\ ~ — ^2 ■ 



V'+,±i/i;ia( x ) = ' i P±v i ;ia( x ) j 



(4.26) 
(4.27) 



eq. ( 4.16| ) becomes: 



£ 



e^^V+.T^l^C^-s j/3 (rr) — (L — > i?) + h.c. (4.28) 



In a similar way the term corresponding to the vector condensate in the la- 
grangian can be written as follows: 



£ 



(V) 



aw 



al^^U^Civ ■ n)^, ]P {x) -(L^R) + h.c. (4.29) 



where n = q/\q\ is the direction corresponding to the total momentum carried 
by the Cooper pair and we have used i/i^Ca ■ nifj + = v ■ ni$LCty+ . 
In the usual basis of the HDET the effective lagrangian is 



£ 



£o + £^ + 



£ 5 A b V ■ d 



A 



A as £ F • d 



= EE*'" 

Here the matrix Aab is as follows: 

A AB = (AorB = 4or5) 
while for A, B = 0, 3 we have: 



t 

AB 



with 



B 



Ar 



( A 













-A 


-£Ai 








+£Ai - 


-A 





V Ai 








-A / 




Ai 


= V 


nA^ . 



X 



B 



(4.30) 



(4.31) 



(4.32) 



(4.33) 



Let us now discuss the precise meaning of the average over velocities. As 
discussed above we embodied a factor of 1/2 in the average over velocities; 
however in the LOFF case the sum over the velocities is no longer symmetric 
and now reads: 
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if k r is kinematical factor of the order of 1 induced by the approximation of the 
Riemann-Lebesgue lemma |34|] and 

(4.35) 



v = (sin^cos^), sin^gsin^), cos9 q ) , 
and 9 q given in ( |4.23| ). For future reference we observe that 

5aS 



^ ViVj = k R 



( sin 2 9„ rrx cos 1 



V 4 



'J3 



(4.36) 



The effective action for the fermi fields in momentum space reads: 



S 



5 f rl A P' 



(4.37) 



u A,S=0 

where (£',£) is the inverse propagator, given by: 



D- A l B {Z',l) 



V ■ £5 AB A AB 



5 4 {£' - £) 



Aab V- 

From these equations one can derive the quark propagator, given by 

/ V-£5 AC A AC \ 

D AB (£j") = (27r)^ 4 (£-£")><Yl 



(4.38) 



c 



Dcb(£) D CB (£) 

A AC V-£6ac 
\ b CB {£) d cb (£) J 



(4.39) 



where 



D CB (£) = (V -IV- 1 - AA f ) 
D CB {£) = (v -£V ■£ - A^A] 



J CB 
)c.B 



(4.40) 



The propagator for the fields x 4 ' 5 does not contain gap mass terms and i 
given by 

{v-iy 1 o 



IS 



D(£,£') = (2vr) 4 5 A {£ - £') 



(V-£) 



-i 



(4.41) 



For the other fields x A i A = 0, • • • ,3, it is useful to go to a representation 
where AA^ and A^A are diagonal. It is accomplished by performing a unitary 
transformation which transforms the basis \ A t ne new basis \ A defined by 



x A = Rabx b , 



with 



RAB = T2 



/ 1 1 \ 

1 -i 

+i -1 

V 1 -1 J 



(4.42) 



(4.43) 



6G 



In the new basis we have 



(AAt 
(At A 



AB 
AB 



where A^, A^ can be approximated as 

X A = ((A + Ai) 2 , (A - Ai) 2 , 



^A^AB 

Xa^ab 



(A + Ai) 2 , (A 



(4.44) 



At) 2 ! 



\a = ((A - AO 2 , (A - Ai) 2 , (A + Ai) 2 , (A + Ax) 2 ) 
For further reference we also define 

He = (A + Ai, Ai - A , A + Ai, Ai - A ) . 



As shown in [34] the integration measure for the LOFF case is 



dH 



2vr / u 2 



(2vr)4 (2tt)3 J_ s "I 



din 



oo 



°°d£o 
2vr 



(4.45) 



(4.46) 



(4.47) 



Let us conclude this section by deriving the gap equation in the LOFF phase. 
As in the CFL and 2SC cases one writes a truncated Schwinger-Dyson equation. 
One obtains two equations: 



A 
Ai 



^ 2 k R 
24A 2 vr 3 



A=0 



24A 2 vr 3 



where 



1(A) 



d 2 £ 



V -IV -£-A 2 + 



i 2ir arcsinh — - 
A 



The values of the parameters have been fixed in paragraph 
equations ( 2.183| ) have a non trivial solution; one gets [34]: 



A f 



Ai 



sinh 



/ 3A 2 7T 2 



cos e q A^ = 



(4.48) 

(4.49) 
The coupled 

(4.50) 
(4.51) 



These results are compatible with the findings of |35] where the scalar conden- 
sate is in the range of a few MeV and A^ is found to be negligible. Since in 
HDET one has Ai = it can be argued that the result Ai ^ found in |35| 

1 



comes from a non leading O 



/'• 



contribution. 
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4.4 Phonon-quark coupling 



Using the results of the two previous sections one gets the following 3-point and 
4-point phonon-quark-quark couplings: 

k 



C 



en + v ■ n 



nA^a. 



and 



V 

if 

-(L^R) + /i.e. (4.52) 
^ ^ [A(% + u • nA^] ^wC^-e 



l ± Y^j}. [y^-n(n- V0)] AW^e^^^C^,^ 



4/V 

xA^aj^^C il> jA _t - (L - i?) + /i.e. 
In the basis of the x fields one gets 



ffl + g\ \~b 
93+94 



(4.53) 



(4.54) 



Here 



93 

94 



icj)A 



AB 



f 

fA AB 
2/2 



+ a AB d[4>} 



+ a AB [^6[<f>] + Q[<t>] 



(4.55) 



with 



1 



0[<f>] = Wq v 
AW 



V</> - n(n • V0) A^ , 

/ 3(n • V<^>) 2 - |V(^| 2 ) - (v • V<p)(n ■ Vcf>) 



, (4.56) 



(4.57) 



/ -1 \ 

-i 
+t 

V +i o o o / 

Finally the effective action for the field (ft, S[4>], is obtained by the lagrangian 
as follows 

/+oo 
dtdxdy- Y £(^(t,iE,j/,fer/g) . (4.58) 

^ fc=— oo 
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4.5 Dispersion law for the phonon field 



To introduce formally the NGB in the theory one uses the same gradient expan- 
sion discussed in the paragraph 2.3.1. At the lowest order one has to consider 
the diagrams in Fig. 2.1, i.e. the self-energy and the tadpole diagrams. At 
the second order in the momentum expansion one gets, after averaging over the 
Fermi velocities 



n(p) 



(J?k R 



vl 



2/2 , 2\ 22 

V±(P X +Py)-V ll P z 



4vr 2 / 2 

One obtains canonical normalization for the kinetic term provided 

H 2 k R 



.r 



2tt 2 



On the other hand, at the lowest order in A^: 



(4.59) 



(4.60) 



— sin' 
2 



8 a + (1 -3cos 2 o 



cos 2 6. 



q ■ 



In conclusion, the dispersion law for the phonon is 

E(p) = J 'v 2 ± {pl + pi) + v^p% 



(4.61) 
(4.62) 

(4.63) 



which is anisotropic. Besides the anisotropy related to v± ^ v\\, there is another 
source of anisotropy, due to the fact that p z , the component of the momentum 
perpendicular to the planes (4.5), differently from p x and p y is a quasi momen- 
tum and not a real momentum. 
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Chapter 5 

Astrophysical implications 



5.1 Overview 



Should we look for a laboratory to test color superconductivity, we would face 
the problem that in the high energy physics programmes aiming at new states 
of matter, such as the Quark Gluon Plasma, the region of the T — ji plane under 
investigation is that of low density and high temperature. On the contrary we 
need physical situations characterized by low temperature and high densities. 
These conditions are supposed to occur in the inner core of neutron stars, under 
the hypothesis that, at the center of these compact stars, nuclear matter has 
become so dense as to allow the transition to quark matter. A schematic view 
of a neutron star is in Fig. 5.1. In the subsequent sections we shall give a 




.7 10.310.6 



core superfluid n crust 



Fig. 5.1 Schematic view of a neutron star as computed by an equation of state 



with three nucleon interaction [75] 



pedagogical introduction to the physics of compact stars and we shall review 
some possible astrophysical implications of the color superconductivity. 
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5.2 A brief introduction to compact stars 



To begin with, let us show that for a fermion gas high chemical potential means 
high density. To simplify the argument we assume that the fermions are mass- 
less and not interacting, so that p = e F = pp. It follows from eq. ( |1,8| ) that 

pocju 3 , (5.1) 

which means that the chemical potential increases as p 1 ' 3 : This is the reason 
why we should search color superconductivity in media with very high baryonic 
density. 

The equations ( |1.8| )-( |i~Tc| ) allow to determine the equation of state, which 
is obtained by solving the equations 

P = $(x F ), poc4, (5.2) 

where xp = pp/mc. For the simplest case (free massless fermions) we obtain 

P oc p 4 , p oc p 3 , (5-3) 

and therefore the equation of state has the form 

P = Kpl . (5.4) 

Clearly this result also holds for massive fermions in the Ultra Relativistic (UR) 
case where the mass is negligible. One can prove that in the Non Relativistic 
(NR) case one has P = Kps. More generally the equation of state can be 
approximated by the expression 

P = Kpi , (5.5) 
and the two cases discussed above are characterized as follows: 

p <C W^g/cm 3 {electrons} 

p <C 10 15 g/cm 3 {neutrons} 

p 3> 10 6 gj cm , {electrons} 

p 3> W 15 g/cm 3 {neutrons} . 

Let us note explicitly that, at T = 0, P ^ 0. This is a quantum- mechanical 
effect due to the Pauli principle and the Fermi Dirac statistics (for comparison, 
for a classical Maxwell Boltzmann gas P — > when T — > 0). In absence of 
other sources of outward pressure it is the pressure of the degenerate fermion 
gas that balances the gravity and avoids the stellar collapse. 



(-§) 



UR (y - 1) 
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One can see that the densities that can be reached in the compact stars are 
very different depending on the nature of the fermions. The two cases corre- 
spond to two classes of compact stars, the white dwarfs and the neutron stars. 
White dwarfs (w.d.s) are stars that have exhausted nuclear fuel; well known 
examples are Sirius B, or 40 Eri B. In the Hertzsprung-Russel diagram w.d.s 
fill in a narrow corner below the main sequence. In a w.d. stellar equilibrium 
is reached through a compensation between the inward pressure generated by 
gravity and the outward the pressure of degenerate electrons. Typical values of 
the central density, mass and radius for a w.d. are p = 10 6 g/cm 3 , M ~ M Q , 
R ~ 5, 000km. 

Suppose now that in the star higher values of p are reached. If p increases, 
inverse beta decay becomes important: 

e~p — > nv . (5-6) 

This process fixes the chemical composition at equilibrium 

Pe + P P = Pn ■ (5.7) 

On the other hand one has to enforce neutrality: 

n e = n p . (5.8) 

One proves that these two conditions fix the ratio of the number of proton n p 
to the number of neutrons n n for ultrarelativistic particles: 

^L = -. (5.9) 

This number should be seen as a benchmark value, as it is derived under sim- 
plifying hypotheses, most notably the absence of interactions and the neglect 
of masses. In any event it suggests that, for higher densities, the star tends to 
have a relatively larger fraction of neutrons and therefore it is named a neutron 
star. It must be stressed that one of the relevant facts about neutron stars is 
that the general relativity effects cannot be ignored and the relevant equilib- 
rium equations to be used are the Oppenheimer-Volkov equations of hydrostatic 
equilibrium. 

The following simple argument, due to Landau (1932) can be used to eval- 
uate the relevant parameters of white dwarfs and neutron stars (see the text- 
book Jr^l ; more recent reviews of compact stars are in JT6|). Let us consider N 

fermions in a sphere of radius R at T = 0; the number of fermion per volume 

„ 1 
unit scales as n ~ N/R ; the volume per fermion is therefore ~ — and the 

n 

uncertainty on the position is of the order of n" 1 / 3 ; the Fermi momentum is of 
the order of the uncertainty on the fermion momentum and therefore 

p F ~ n 1/3 h , 

a result we obtained already under more stringent hypotheses (Fermi-Dirac 
distribution) and derived again here using only the uncertainty relations. The 
Fermi energy of the baryons is therefore 

hcN 1 / 3 
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if N is the total number of baryons. Note that this applies both to neutron 
stars and to electron stars, because also in stars where the pressure mainly come 
from electrons there will be a considerable amount of protons and neutrons and 
the largest part of the energy comes from the baryons, not from the electrons. 
On the other hand the gravitational energy per baryon is 

GNm 2 B 
R — ' 

and the total energy can be estimated as 

hcN 1 / 3 GNm% 
E = E G + E F jj-B . (5.10) 

Now equilibrium can exist only if E > 0. As a matter of fact if E < (JV large) 
then lirriR^o E = — oo, which means that the energy is unbounded from below 
and the system is unstable. Therefore, E > gives the maximum number of 
baryons as follows: 

/ he \ 3/2 

N < N max = y^r) ~ 2 x 10 57 . (5.11) 

As a consequence the maximum mass is 

M max = N max m B = 1.5M . (5.12) 

This mass can be estimated better and its better determination (~ 1.4M Q ) 
is known as the Chandrasekhar limit; for our purposes the estimate ( 5.12[ ) is 



however sufficient. Notice that the Chandrasekhar limit is similar for compact 
stars where the degeneracy pressure is mainly supplied by electrons and those 
where it is supplied by baryons. 

One can also estimate the radius of a star whose mass is given by ( 5.12| ). 
One has 

e F >mc 2 ~% NUL ~ % ( V2 (5.13) 



R " max R \Gm 2 B 



and, therefore, 

fl f He \ 1/2 r 5x 10 s em {m=m e } 



R 



) -i 



mc \Gm 2 R J I 3xl0 5 cm {m=m n }. 



If a neutron star accretes its mass beyond the Chandrasekhar limit nothing 
can prevent the collapse and it becomes a black hole[|. In the following table 
we summarize our discussion; notice that we report for the various stars also 
the value of the parameter GM/Rc 2 i.e. the ratio of the Schwarzschild radius 
to the star's radius. Its smallness measures the validity of the approximation 
of neglecting the general relativity effects; one can see that for the sun and 
the white dwarfs the newtonian treatment of gravity represents a fairly good 
approximation . 

The exact determination of the mass limit depends on the model for nuclear forces; for 
example in [R7| the neutron star mass limit is increased beyond IAMq. 
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M 


R 


\cm° / 


GM 


Sun 


M & 




1 


10~ 6 


White Dwarf 


<M 




< 10 7 


1(T 4 


Neutron Star 


1 - 3M© 


10 _5 i? o 


< 10 15 


ht 1 


Black Hole 


arbitrary 


2GM 

c 2 


M 


1 



Table 5.1 Paramers of different stellar objects. 



Neutron stars are the most likely candidate for the theoretical description 
of pulsars. Pulsars are rapidly rotating stellar objects, discovered in 1967 by 
Hewish and collaborators and identified as rotating neutron stars by Gold |78tl ; 
so far about 1200 pulsars have been identified. 

Pulsars are characterized by the presence of strong magnetic fields with 
the magnetic and rotational axis misaligned; therefore they continuously emit 
electromagnetic energy (in the form of radio waves) and constitute indeed a 
very efficient mean to convert rotational energy into electromagnetic radiation. 
The rotational energy loss is due to dipole radiation and is therefore given by 



dE 
~dt 



du B 2 R 6 lu 4 sin 2 9 



dt 



6c 3 



(5.14) 



Typical values in this formula are, for the moment of inertia I ~ R 5 p 
10 45 g/cm 3 , magnetic fields B ~ 10 12 G, periods T = 2tt/lo in the range 1.5 

msec-8.5 sec; these periods increase slowly 0, with derivatives 



dT 
~dt 



10 



-12 



10~ 21 , and never decrease except for occasional jumps (called glitches). 

Glitches were first observed in the Crab and Vela pulsars in 1969; the vari- 
ations in the rotational frequency are of the order 10 -8 — 10 -6 . 

This last feature is the most significant phenomenon pointing to neutron 
stars as a model of pulsars in comparison to other form of hadronic matter, 
such as strange quarks. It will be discussed in more detail in par. ^1], where we 
will examine the possible role played by the crystalline superconducting phase. 
In the subsequent three paragraphs we will instead deal with other possible 
astrophysical implications of color superconductivity. 



2 Rotational period and its derivative can be used to estimate the pulsar's age by the 
approximate formula dT , see e.g. [79[ . 
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5.3 Supernovae neutrinos and cooling of neutron stars 



Neutrino diffusion is the single most important mechanism in the cooling of 
young neutron stars, i.e. with an age < 10 5 years; it affects both the early 
stage and the late time evolution of these compact stars. To begin with let us 
consider the early evolution of a Type II Supernova. 

Type II supernovae are supposed to be born by collapse of massive (M ~ 
8 - 2OM ) stars 3 . These massive stars have unstable iron cores0 with masses 
of the order of the Chandrasekhar mass. The explosion producing the super- 
nova originates within the core, while the external mantle of the red giant star 
produces remnants that can be analyzed by different means, optical, radio and 
X rays. These studies agree with the hypothesis of a core explosion. The emit- 
ted energy (~ 10 51 erg) is much less than the total gravitational energy of the 
star, which confirms that the remnants are produced by the outer envelope of 
the massive star; the bulk of the gravitational energy, of the order of 10 53 erg, 
becomes internal energy of the proto neutron star (PNS). The suggestion that 
neutron stars may be formed in supernovae explosions was advanced in 1934 



by Baade and Zwicky [SO] and it has been subsequently confirmed by the ob- 
servation of the Crab pulsar in the remnant of the Crab supernova observed in 
China in 1054 A.D. 

We do not proceed in this description as it is beyond the scope of this review 
and we concentrate our attention on the cooling of the PNS f], which is mostly 
realized through neutrino diffusion. By this mechanism one passes from the 
initial temperature T ~ 20 — 30 MeV to the cooler temperatures of the neutron 
star at subsequent stages. This phase of fast cooling lasts 10-20 sees and the 
neutrinos emitted during it have mean energy ~ 20 MeV. These properties, 
that can be predicted theoretically, are also confirmed by data from SN 1987A. 

The role of quark color superconductivity at this stage of the evolution of 
the neutron stars has been discussed in fl82| . In this paper the neutrino mean 
free path is computed in a color super-conducting medium made up by quarks in 
two flavor (2SC model). The results obtained indicate that the cooling process 
by neutrino emission slows down when the quark matter undergoes the phase 
transition to the superconducting phase at the critical temperature T c , but then 
accelerates when T decreases below T c . There should be therefore changes in 
the neutrino emission by the PNS and they might be observed in some future 
supernova event; this would produce an interesting test for the existence of a 
color superconducting phase in compact stars. 

Let us now consider the subsequent evolution of the neutron star, which also 
depends on neutrino diffusion. The simplest processes of neutrino production 
are the so called direct Urea processes 



3 The other supernovae, i.e. type I supernovae, result from the complete explosion of a star 
with 4M Q < M < SM Q with no remnants. 

4 Fusion processes favor the formation of iron, as the binding energy per nucleon in nuclei 
has a maximum for A ~ 60. 

5 See [M for further discussions. 
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/2-/i+*-/2 + ^; (5.15) 

by these reactions, in absence of quark superconductivity, the interior temper- 
ature T of the star drops below 10 9 K (~ 100 KeV) in a few minutes and in 
10 2 years to temperatures ~ 10 7 K. Generally speaking the effect of the for- 
mation of gaps is to slow down the cooling, as it reduces both the emissivity 
and the specific heat. However not only quarks, if present in the neutron star, 
but also other fermions, such as neutrons, protons or hyperons have gaps, as 
the formation of fermion pairs is unavoidable if there is an attractive attraction 



in any channel (see the discussion in par. L3). Therefore, besides quark color 
superconductivity, one has also the phenomenon of baryon superconductivity 
and neutron superfluidity, which is the form assumed by this phenomenon for 
neutral particles. The analysis is therefore rather complicated; the thermal 
evolution of a late time neutron star has been discussed in [ p3|j , but no clear 
signature for the presence of color superconductivity seems to emerge from the 
theoretical simulations and, therefore, one may tentatively conclude that the 
late time evolution of the neutron stars does not offer a good laboratory to test 
the existence of color superconductivity in compact stars. 



5.4 i?-mode instabilities in neutron stars and strange 
stars 

Rotating relativistic stars are in general unstable against the rotational mode 
(r-mode instability) |84[| . The instability is due to the emission of angular 
momentum by gravitational waves from the mode. Unless it is damped by 
viscosity effects, this instability would spin down the star in relatively short 
times. More recently it has been realized [^] that in neutron stars there is 
an important viscous interaction damping the r-mode i.e. that between the 
external metallic crust and the neutron superfluid. The consequence of this 
damping is that r-modes are significant only for young neutron stars, with 
periods T < 2 msec. For larger rotating periods the damping of the r— mode 
implies that the stars slow down only due to magnetic dipole braking. 

All this discussion is relevant for the nature of pulsars: Are they neutron 
stars or strange stars? 

The existence of strange stars, i.e. compact stars made of quarks u, d, s in 
equal ratios would be a consequence of the existence of stable strangelets. This 
hypothetical form of nuclear matter, made of a large number of u, d, s quarks 
has been suggested by Bodmer and Witten [^] as energetically favored in 
comparison to other hadronic phases when a large baryonic number is involved. 
The reason is that in this way the fermions, being of different flavors, could 
circumvent the Pauli principle and have a lower energy, in spite of the larger 
strange quark mass. If strangelets do exist, basically all the pulsars should 
be strange stars because the annihilation of strange stars, for example from a 
binary system, would fill the space around with strangelets that, in turn, would 
convert ordinary nuclear stars into strange stars. 
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An argument in favor of the identification of pulsars with strange stars is 
the scarcity of pulsars with very high frequency (T < 2.5 msec). 

This seems to indicate that indeed the r— mode instability is effective in slow- 
ing down the compact stars and favors strange stars, where, differently from 
neutron stars, the crust can be absent. Even in the presence of the external 
crust, that in a quark star can be formed by the gas after the supernova explo- 
sion or subsequent accretion, the dampening of the r — mode is less efficient. As 
a matter of fact, since electrons are only slightly bounded, in comparison with 
quarks that are confined, they tend to form an atmosphere having a thickness 
of a few hundred Fermi; this atmosphere produces a separation between the 
nuclear crust and the inner quark matter and therefore the viscosity is much 
smaller. 

In quark matter with color superconductivity the presence of gaps A S> T 
exponentially reduces the bulk and shear viscosity, which renders the r—mode 
unstable. According to Madsen [88, (8iJ this would rule out compact stars 
entirely made of quarks in the CFL case (the 2SC model would be marginally 
compatible, as there are ungapped quarks in this case). For example for A > 1 
MeV any star having T < 10 msec would be unstable, which would contradict 
the observed existence of pulsars with time period less than 10 msec. 

However this conclusion does not rule out the possibility of neutron stars 
with a quark core in the color superconducting state, because as we have stressed 
already, for them the dampening of the r—mode instability would be provided 
by the viscous interaction between the nuclear crust and the neutron superfluid. 



5.5 Miscellaneous results 

Color superconductivity is a Fermi surface phenomenon and as such, it does 
not affect significantly the equation of state of the compact star. Effects of 
this phase could be seen in other astrophysical contexts, such as those con- 
sidered in the two previous paragraphs or in relation to the pulsar glitches, 
which will be examined in the next paragraph. A few other investigations have 
been performed in the quest of possible astrophysical signatures of color super- 
conductivity; for instance in (5^] it has been suggested that the existence of a 
2SC phase might be partly responsible of the gamma ray bursts, due to the 
presence in the two-flavor superconducting phase of a light glueball that can 
decay into two photons. Another interesting possibility is related to the stabil- 
ity of strangelets, because, as observed in J90"[ , CFL strangelets, i.e. lumps of 
strange quark matter in the CFL phase may be significantly more stable than 
strangelets without color superconductivity. 



Finally I wish to mention the observation of [ 11 ] concerning the evolution of 
the magnetic field in the interior of neutron stars. Inside an ordinary neutron 
star, neutron pairs are responsible for superfluidity, while proton pairs produce 
BCS superconductivity. In this condition magnetic fields experience the ordi- 
nary Meissner effect and are either expelled or restricted to flux tubes where 
there is no pairing. In the CFL (and also 2SC) case, as we know from paragraph 
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3.1.2| , a particular U(l) group generated by 



Q = 1®Q + T®1 (5.16) 

remains unbroken and plays the role of electromagnetism. Instead of being 
totally dragged out or confined in flux tubes, the magnetic field will partly 
experience Meissner effect (the component Aa), while the remaining part will 
remain free in the star. During the slowing down this component of the magnetic 
field should not decay because, even though the color superconductor is not a 
BCS conductor for the group generated by (5.16), it may be a good conductor 



due to the presence of the electrons in the compact star. Therefore it has been 
suggested | ]] that a quark matter core inside a neutron star may serve as an 
"anchor" for the magnetic field. 



5.6 Glitches in neutron stars 

Glitches are a typical phenomenon of of the pulsars, in the sense that probably 
all the pulsar have glitches (for a recent review see pl| ). Several models have 
been proposed to explain the glitches. Their most popular explanation is based 
on the idea that these sudden jumps of the rotational frequency are due to the 
the angular momentum stored in the superfluid neutrons in the inner crust (see 
Fig. 5.1), more precisely in vortices pinned to nuclei. When the star slows 
down, the superfluid neutrons do not participate in the movement, until the 
state becomes unstable and there is a release of angular momentum to the 
crust, which is seen as a jump in the rotational frequency. 

The presence of glitches is one of the main reasons for the identification 
of pulsars with neutron stars; as a matter of fact neutron stars are supposed 
to have a dense metallic crust, differently from quark stars where the crust is 
absent or, if present, is much less dense (~ 10 11 g cm" 3 ). 

A schematic and a somehow exaggerated diagram of glitches is shown in 
fig. 5.2. The interesting aspect of the LOFF phase discussed in chapter || is 
that even in quark stars, provided we are in a color superconductivity phase, 
one can have a crystal structure able to produce glitches; it would be given 
by a lattice characterized by a geometric array where the gap parameter varies 



periodically. This would avoid the objection raised in 92] that excludes the 



existence of strange stars because, if the strange matter exists, strange stars 



should be rather common, as we discussed in par. 5.4, in contrast with the 
widespread appearance of glitches in pulsars. Therefore, if the color crystalline 
structure is able to produce glitches, the argument in favor of the existence of 
strange stars would be reinforced. 

In a more conservative vein one can also imagine that the LOFF phase 
be realized in the inner core of neutron star; in this case the crystalline color 
superconductivity could be partly responsible for the glitches of the pulsar. 
Some estimates of the effect are discussed in ]35]l ; a detailed analysis of this 
scenario is however premature as one should first complete the study of the 
LOFF phase in two directions, first by including the third quark (all the study 
in section |] was relative to the two flavor case) and, second, by sorting out the 
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Fig. 5.2 Schematic view of the period P of the PSR 0833-45 (Vela) showing 
four glitches in the years 1969-1980. Other glitches have been observed in the 
subsequent years. 

exact form of the color lattice, as the plane wave ansatz discussed in section |I| 
might not correspond to the true vacuum of the theory. 
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